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PREFACE 


This book is an attempt to present the field of analogue computa¬ 
tion and simulation in a compact and handy form. There is a 
large number of books, almost all American, which present 
detailed, thorough, intricate, and all-embracing courses on ana¬ 
logue computing, with the aid of which the student could become 
a professor of computation, mathematics, control theory and 
electronics all rolled into one. I have started with the assumption 
that the reader is intelligent enough to realize the extent and scope 
of a method without being given large numbers of detailed and 
unreadable examples; and that, having understood the elements 
of the subject from the bookwork and examples given, he can 
use his ability to read specialized works and to chase references. 

The reader is assumed to have a knowledge of the calculus, 
mathematics and physics such as would be adequate to gain ad¬ 
mission to a university. A certain amount of mathematics is 
presented in the first chapter, omission of which will not greatly 
impair understanding of the remainder of the work. However, in 
view of the widespread confusion I have found to exist on the 
topics involved, it was thought well worth including this back¬ 
ground. 

The electronic equipment of an analogue computer is described 
in Chapter 2, and its use to solve simple problems shown in 
Chapter 3, together with the method of scaling. Chapter 4 
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describes the general layout of a computer with a description of 
some commercial computers. The use of a large computer is 
described in Chapter 5 by means of a particular example—a 
reactor simulation. Types of problem which have been simulated 
on large computers are described briefly. Finally, Chapter 6 out¬ 
lines some general analogue devices. The differential analyser is 
described for historical reasons and because of the insight it has 
been found to give some people into the solution of equations by 
simulation methods. The methods used in the incremental com¬ 
puter (or digital differential analyser) are explained briefly. Con¬ 
ductive analogues and some other methods are described, and, 
finally, some special techniques, which have been used in con¬ 
nexion with electronic analogue computing, are given. 

I am very greatly indebted to the authors of many works, 
notably W. J. Karplus’s excellent book Analogue Simulation in 
the Solution of Field Problems, and also to A. E. Rogers and 
T. W. Connoly, Analog Computation in Engineering Design, both 
published by McGraw-Hill. Analog Computation, by A. C. 
Jackson (McGraw-Hill), is an extremely useful book. The biblio¬ 
graphy included at the end of the book will be helpful to those 
who would dig deeper into the subject than this book attempts. 

I would like to express my gratitude to Percy Hammond for 
asking me to write this book, and to Dr. T. O. Jeffries, under 
whom I was working at the time, who gave me much useful 
advice and encouragement with the text. 

My thanks are due to Solartron and to Pace systems for per¬ 
mission to publish description of their computers, and to the 
English Electric Co. for permission to write this book. 

Finally, I would like to record my thanks to Miss I. London 
who typed the manuscript. 




chapter 1 

INTRODUCTION 

AND 

MATHEMATICAL 

TECHNIQUES 


INTRODUCTION 

Geometry is one of the oldest branches of mathematics. It is 
supposed to have arisen out of the desire of the Ancient Egyptians 
to make accurate maps in order to settle disputes over land¬ 
holders’ boundaries. Their maps were analogues of the spatial 
relationships of neighbouring areas of land, and the analogue was 
available in a convenient form after the floods of the Nile had 
removed the boundary markers. 

Today every draughtsman is familiar with the geometrical 
methods first thought out by the Greeks and Egyptians. With their 
aid he can lay out in a convenient and readily changed form the 
design of a complex machine or building to high accuracy. 

Man does much of his thinking by drawing analogies between 
something which is partly understood and something which is 
more fully understood or which is more readily handled. The 
physicist explains some aspects of the behaviour of atoms and 
electrons by comparing them to billiard balls, the barrister points 
to a precedent, a case similar to the one in question, to help decide 
the answer to his case, and the theatre seeks to clarify our personal 
relationships by acting a similitude of real life. 

Many devices have been thought up to help in the understand¬ 
ing of the problems of engineering, physics and mathematics. 
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Mathematics and physics are concerned, among other things, with 
finding models, expressed as relationships between numbers, for 
the world as our senses perceive it. And the engineer applies these 
models, which are analogues of the physical world. 

This book is intended to give an outline of some of the many 
models and analogues that have been produced to solve physical 
problems for the engineer, and, more particularly, to show how to 
use and programme the modern electronic analogue computer. 


ANALOGUE AND DIGITAL COMPUTING 
MACHINES 

Basically the digital computer does sums, sums that most 
people could do by the age of 10 or so, takes two or three simple 
types of decision, has an enormous memory, and works at very 
high speed. It does arithmetic. 

The analogue computer is any arrangement which can be made 
to behave in as near the same way as possible to some real system 
of some kind. Most, if not all, digital computers are essentially 
serial in operation (even those known as parallel machines), that 
is, the numbers stored in the machine appear in some central 
arithmetic unit one after another and are then added to or com¬ 
pared with some other number before being sent back to the 
memory store. At one time, only one “sum” is being done on the 
whole set of numbers in the computation. By the use of very small 
computing times, of the order of a few microseconds, and stores 
of very large capacity, very large and complex problems can be 
handled. 

Analogue machines are essentially parallel devices. All the 
numbers describing the problem, represented as voltages or shaft 
positions or as pressures, flows or distances, are all processed 
simultaneously. Additions, multiplications and integrations are 
all proceeding at the same time. A widely used device is the wind 
tunnel where the physical variables such as air speed and force are 
represented by air speed and force, but are scaled down suitably to 
give a representation of a large structure by a small model. 

The digital machine is capable of great accuracy and speed, but 
it needs a skilled mathematician to programme it. Great care must 
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be taken to avoid giving the machine impossible situations, or 
getting it going round in circles trying to find an answer to a 
faulty logical situation. On an analogue, an engineer can quickly 
get the feel of a design and alter parameters and see what happens. 
He will often operate the machine himself. 

The experts say that any problem can be solved on a digital 
machine. But it may take far longer to programme, test the pro¬ 
gramme, and run off the solutions than to run an analogue 
machine programme. Nuclear reactor studies, rocket systems, and 
problems involving many non-linearities can frequently be solved 
most economically on an analogue computer. It is possible to 
include parts of the physical system under test in the analogue 
set-up in some cases. 

The accuracy of an analogue machine is seldom better than one 
part in 1000. This is better than the accuracy of the physical data 
of many problems. But this is sometimes too low, and the 
additional cost of a digital programme becomes necessary. A 
fairly detailed costing of the computation of some radiation 
integrals involving Bessel functions has been done by Hansen. 
The analogue solution, accurate to two figures, took 2 hours to 
programme, 50 minutes to compute the integrals required, and 
was costed at $53. The digital programme took 2 weeks to write 
and 50 minutes on the machine, and was costed at $1377. This 
was much the same as the cost of a hand calculation. 

However, at the present time, a large analogue problem can tie 
up a computer for weeks and months, and no other problem can 
be put on the machine during this time since the constants of the 
problem are represented as potentiometer settings. This can be 
overcome if some automatic means is used for setting up the 
problem, but these are not yet widely available. In contrast, a 
digital computer reads in its programme and data, computes, and 
if everything goes well, prints out the answers and is clear for the 
next operator in 10 minutes, an hour, or in extreme cases a night’s 
work. 


HISTORICAL 

The first device described for solving differential equations was 
an integrating machine developed from the principle of the plani- 
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meter by Prof. James Thomson, Lord Kelvin’s brother. Its con¬ 
struction and use was described in a paper to the Royal Society 
in 1875. A ball is allowed to roll on a disc, and the number of 
rotations is transmitted to a roller. This is described in Chapter 5. 
The roller then gives a total rotation proportional to the integral 
of the distance of the ball from the centre of the disc with respect 
to the rotation of the disc. 

Although Lord Kelvin indicates how the device can be used for 
the solution of differential equations, and quotes problems of the 
vibration of stretched chords, water waves in channels, heat trans¬ 
fer and the motion of the tides as being suitable, it was not used 
until Vannevar Bush re-invented an integrator at M.I.T. in 
America, and built a six-integrator machine in the 1920’s. This 
and a later machine were used in the solution of problems of 
synchronous motor control and stability, and emitting filaments, 
among many others. During and after the Second World War, 
various machines were built, ranging from a two-integrator 
Meccano machine built by Prof. Hartree, accurate to l\ per cent, 
to a thirty-integrator machine, accurate to 1 part in 10,000, built 
at M.I.T. with electric couplings and drive, and a paper tape 
input and output. These machines were widely used for such 
things as control and vibration problems. 

After the Second World War, valves were good enough and 
cheap enough for the construction of high-gain d.c. amplifiers of 
fairly good stability to be possible. Several computers were built 
using these amplifiers to perform summing and integrating 
functions in a manner that will be described later. Careful design 
enabled reliable computers of up to about 100 amplifiers to be 
built, but the difficulty of amplifier output drift was only overcome 
by the widespread introduction of the chopper relay stabilization 
technique developed by E. A. Goldberg in about 1950. Since 
then, a great many computers of up to about 100 amplifiers and 
some of as many as 1500 amplifiers have been built. These are 
reliable, accurate to about 0-1 per cent, and have been used for 
an enormous variety of problems in, for example, control, heat 
transfer, chemical engineering, missile system development, and 
nuclear reactor studies. 

The development of analogue computers has been parallel to 
and in some ways complementary to the development of digital 
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computers. A great deal of nonsense and rivalry has appeared, 
where one method has been acclaimed as so superior to the other 
that the other is considered to be outdated. But each method of 
approach to the solution of a problem will have its uses. Analogue 
solutions are often cheaper, quicker, simpler to programme, but 
far less accurate than digital computer solutions; but there are 
many problems which it is not reasonable to attempt to put on to 
an analogue machine, such as the solution of large sets of linear 
equations, the accurate evaluation of integrals, or the compilation 
of tables, e.g. logarithms. The fact that both methods have their 
own uses is shown by the many firms which have found it economic 
to install both analogue and digital computing systems. Often one 
will feed data to the other, and particular cases can be solved on 
both, enabling cross-checking of results to be carried out. 


MATHEMATICAL TECHNIQUES 

The programmer of an analogue computer may be given sets 
of equations to solve which involve unfamiliar notations and 
mathematical language, although the equations themselves do not 
involve any really difficult mathematics. The writer has met many 
users of analogue machines who have been puzzled by the mathe¬ 
matical language of the problem which they are in fact solving. 
Therefore, the following is provided as a sort of mathematical 
“phrase book”, without any attempt at a serious study of the 
“grammar” of the mathematics involved. It is hoped that the 
trained mathematicians will forgive me for the sketchy way the 
subjects are covered. I know that the material of this chapter 
has proved helpful to my engineering colleagues, but I fear the 
wrath of my mathematical lecturers and teachers. 

On this basis, a short outline of the Laplace transform is given 
and of the use of the operator p (or s). The operator del = V is 
explained briefly, and types of equation where it appears are 
mentioned. A short example illustrates the method of separation 
of variables. Finite difference forms of differentials and Lagrange 
coefficients are given. The notation of matrices is explained, and 
the existence of eigenvalues mentioned, together with some re¬ 
marks on eigenfunctions. 
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THE LAPLACE TRANSFORM 

The Laplace transformation provides a technique for solving 
linear differential equations by algebraic methods. A short out¬ 
line is given here. The student will find references to fuller 
treatments at the end of the chapter. 

We define the transform J(p) of a function/(r) by the relation 

/(P) = j o 7(0 e ~*df (1) 

This is subject to mathematical conditions concerning the limits 
and validity of the integration, which are satisfied in the majority of 
cases ordinarily encountered. Two theorems are readily derived: 

If/O) is the Laplace transform of f(t), 

(1) J(p+a) is the transform of /(t)e" a( ; 

(2) (d/dt)/(r) transforms to pj(p)-y 0 , where y 0 = /(0+),f and 
applying this result twice, 

(d 2 /d t 2 )f(t) transforms to p 1 J(p)-py 0 -y u 


Table 1 


Function 


Laplace transform 

Impulse 

<5(0 

1 

Step 

Hit) 

HP 

Ramp 

t 

l Ip 2 



i 


n 

p n+ 1 


e- 

l 

P + a 


i n e~ at 

(P + a )» 


sin cot 

co 2 +p 2 


d/ 

d t 

pf(p) -/(0 +) 


t The notation /(0 +) denotes the limiting value of /(?) as 1 tends to 
zero, where t is always greater than zero, 
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where and y 0 are the values of df/dt and / at time / = 0 +. 
It can easily be shown that the transform of a sum of two 
functions is the sum of the transforms, and that the transform 
of Kf(t) is Kj(p), where A - is a constant. 

Integration of equation (1) leads to the transforms shown in 
Table 1. 


Exercise 

Prove the above theorems by integration from the definition of the 
transformation. 


Solution of a Simple Equation 

Consider the equation ax+bx+cx = Take the Laplace 
transform of both sides; 

Then ( ap 2 + bp+c)x = kf(p)+(ap+b)x 0 + ax l 


thus 


-z n _kf(p)+(ap+b)x 0 + ax l 
X(P) ap 2 + bp+c 


By expanding in partial fractions and the use of the table of 
transforms, the function x(t) can be found given the initial 
conditions and the function kf(t). This is a simple matter for this 
equation, and the solution can be found as quickly by other means. 
But the technique is as easy to apply for simultaneous linear 
differential equations, or equations of high order, the difficulty 
lying in the resolution into partial fractions of the expression 
equivalent to the r.h.s. in the equation above. 

Some people use the symbol s for the Laplace operator rather 
than p. The writer can detect no majority view, and no reason 
to be partisan. Since the symbols are interchangeable, p has been 
selected on arbitrary grounds. 


THE TRANSFER FUNCTION 

If a system has input /(f), output x(t), then the function 


’ f(p) 


is called the transfer function of the system. Note that it can also 
be obtained by replacing the differential operator d/d t by the 
symbol p in the differential equation of a linear system. 
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Sinusoidal Inputs 

For a linear system with sinusoidal input, all the system vari¬ 
ables will be sinusoidal, differing only in amplitude and phase 
from the input. 

If F(t) = Ae +jo>t = (a + jb) e +J “', 

where a and b are real numbers, then if 

A = a 2 + b 2 and tan (j> = b/a, 
it is easily shown that 

F(t) = A cos (cot + (j>) +jA sin (cot + </>). 

So we can represent a sinusoid in phase and magnitude, in a linear 
system, by a complex number of the form Ae Jat . Any algebra 
involved in the system can be performed on this number. Taking 
real parts at any stage gives the same result as would have been 
obtained if the solution had been worked through by operating 
on the real part alone. Now 

^-(Ae +Ja> ‘) = jo)Ae +Ja>t . 
at 

Therefore the operation of differentiation is formally identical to 
multiplication by yen. It is convenient to write jeo = p, so 
d/d t =p= jco. 

The differential equation describing a linear system can there¬ 
fore be written in the form 

a n P n X+ a„- 1 p" 1 X +... + ai pX+a 0 X 

= b m p m F+b m - l p m ~ 1 F + ... + b l pF+b 0 F, 

where F is the input and X the output sinusoid of the system. 

The transfer function defined as the ratio of the Laplace trans¬ 
forms of output and input would give exactly the same expression 
except that the complex numbers X and F representing the output 
and input functions x and / would be replaced by the transforms 
x and /. The term “transfer function” is used to cover both 
expressions. Care should be taken in interpreting this function. 
When it is used in the sinusoidal sense, the modulus of G(jco) 
will give the gain, and the ratio of its imaginary to real parts the 
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tangent of the phase. When used in its Laplace transform sense, 
it is an operator, operating on the transform of the input to give 
the transform of the output. 

However, since the two expressions are formally identical, we 
can obtain the gain and phase of the transfer function simply by 
replacing the symbol p by jco and taking the modulus and argu¬ 
ment as before. 


Decibels 

The decibel is the measure of the logarithm of the ratio between 
two quantities which measure powers. The relationship in deci¬ 
bels (abbreviated to dB) between powers P 0 and P t is defined by 

number of decibels = 101og 10 

Thus, if the input and output voltages of a system are V i and V 0 , 
since the power is proportional to the voltage squared, the gain 
in decibels of the system is 


gain in decibels = 20log, 0 


Now, consider the simple lag, with transfer function 
G(p)=l/(l+pT). 

The gain is the modulus of this, with p replaced by joo; 

?( “ ) “ |t+J3t| “ 


n decibels = 20log 


V V 1 +*> 2 T 2 J 

ioiog(i+co 2 r 2 ). 


When co is small, the gain is = 0 dB. 

When co is large, the gain is = — 101og(co 2 r 2 ), 

= —201ogo)T. 

Thus for large values of co the gain decreases by 20 dB for each 
decade (factor of 10) in frequency; this is equivalent to 6 dB for 
each octave (doubling in frequency), since log2 = 0-3 very nearly. 
Similarly, the phase is given by = —tan _1 a>T. Now, on 
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a logarithmic scale of m, the points corresponding to to = a/7 
and co-l/aT are symmetrical about the point to= 1 /7; also, 
since at the point to = 1/7, </> = — 45°, and since 

. a ,1 

tan -1 —. 7+ tan _1 — .7 = 90°, 

7 a7 

the curve of phase against frequency, with frequency on a logar¬ 
ithmic scale, is symmetrical about 45° and to = 1/7, as shown in 
Fig. 1. 


Exercises 

1. Calculate the gain and phase of the transfer functions G(p) = p, 
lip, 1 — pT, and 1/(1 + pT) as a function of frequency. Show that 
the gain of the third and fourth can be approximated by appropriately 
placed asymptotes with slope 6dB/octave as shown in Fig. la. 

2. A network has transfer function 

G(P ^ (1 +/>)0 +01/))' 

What is its response to a ramp input ? By using the asymptote method 
outlined in Ex. 1, sketch its frequency response on a log-linear scale. 


The p-plane 

We may factorize a transfer function G(p) to give an expression 
of the form 


G(p) 


(fix -pWi-P)~. 

(ai—p)(a 2 -p)...’ 


where the a’s and P's are real, or complex pairs, or may be re¬ 
peated. 

Then, if p = a + jco, where a and to are real, we may mark on 
the Argand diagram of the complex numbers the poles of G(p), 
where p = & u a 2 ,..., and the zeros, where p = p lt p 2 ,... (A pole 
of G(p) is a point where its value is indeterminately large.) 
Suppose we have a pair of poles at a+jb, a—jb; in the partial 
fraction expansion of the response there will therefore be a term 
of the form 

A B 


p — (a +jb) p — (a —jb )' 
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Reference to the table of transforms shows that the output con¬ 
tains terms of the form Ae a '.e jb '+Be al .e~ Jbt . Thus, if a is 
greater than 0, this represents an unstable response. So poles 
in the right-hand half-plane mean a diverging response. 

Details of the techniques developed for the investigation of 
transfer functions are given in any textbook on control theory. 


Exercise 

Investigate the behaviour of the system with transfer function 


G(p) = 


p 2 + 2£ (On p + (On 2 


as the damping coefficient { varies from + a> to — oo. A thorough 
study of this function will give a great deal of insight into the behaviour 
of many physical systems. 


electrical impedance 

Ohm’s law is well known in the form E = IR. If in the differ¬ 
ential equations for the flow of current through a capacitor and 
an inductance, p is written for d/d t, we obtain E = (\/Cp).I and 
E = Lp.I. The quantities \jCp and Lp are called the impedances 
of a capacitor and an inductance respectively in operational 
form. Reference to the differential equations quickly show that 
the rules for combining impedances in series and parallel are the 
same as those for combining resistors. 


Exercise 

What is the impedance of an inductance, a capacitor and a resistance 
in series, and in parallel? Show that the impedances and transfer 
functions of the networks given in Fig. 2 are those given in the table. 
Derive other transfer functions from combinations of various im¬ 
pedances. 

Kirchhoff’s Laws 

(1) At any junction of impedances, called a node, the total 
inflow of current equals the total outflow. 

(2) Round any closed path in a circuit, called a loop, the sum 
of the currents in each element times its impedance is equal 
to the voltage drop round the loop. 

These follow from the conservation of electric charge at the 



(*>) 



pc 

pL 


+ pCR 
PC 


-Tf> 

R n c 

—vw—1[— 

R, 

|-^VNAA-j 

“Unfr 


~“G3~ '’( l+ 'c7 ! n: ! *) 

Network Transfer Function 


*i(I + P CR 2 ) 

1 + pC(Ri + R 2 ) 


(1 + pC 2 R)(Ci + Ci) 


pCR 
1 + pCR 


, R\Ri 

1 Rl + *2 


Fig. 2. 
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node, and the application of Ohm’s law to each element of the 
loop. 


VECTORS AND SCALARS 

It will be recalled that many physical variables have a direction 
as well as a magnitude associated with them. Force and velocity 
are familiar examples of vectors. Vectors are combined by the 
parallelogram rule, normally introduced as the parallelogram of 
forces. A vector is distinguished in print by bold lettering. 

A scalar quantity is one such as mass or volume which com¬ 
bines by simple addition. Two pounds plus three pounds makes 
five pounds. This is contrasted with vectors, which combine thus: 
a ship heading 30° west of north at 10 knots in a river flowing 
due east at 5 knots will move due north at 5,^/3 knots. 

A vector can be expressed as a function of position; the force 
due to gravity near the earth is directed towards the earth and is 
inversely proportional to the square of the distance from the 
earth’s centre. Similarly, a scalar quantity can be expressed as a 
function of position—for example, the temperature of the air in 
a room as a function of the position in the room. 

The mathematical use of vectors has enabled many branches of 
physics to be described in a compact notation. Differentiation of 
a vector takes a special form, and various related operations of 
differentiation are found convenient. These relate to the import¬ 
ant differential equations such as Laplace’s equation, the diffusion 
equation and the wave equation. 


The Operator V (Del) 

“Del” is a differential operator defined by 


V 


.a .a , a 

1—+ J^- + k r- 

ox oy dz 


where i, j, k, are unit vectors in the x-, y-, and z-directions. 

If <j)(x,y,z) is a scalar function defined in a coordinate system 
( x,y,z ), then the vector 


grad <f> = V0 = i^r+ijr+kjr 
dx 8y dz 


is called the gradient of 4>. It is analogous to the slope of a 
function f(x) that is, f'(x) = (djdx)f(x). The three components 
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of the slope are given by the derivatives in the directions of the 
three axes. 

If a(x,y,z) = (a x ,a y ,a z ) is a vector defined in the system of 
coordinates, then the scalar 


diva = V. 


8a x 8a y 8a z 
ox 8y 8z 


( l l +i | +k l) <i “- +i “’ +k “ J 


is called the divergence of the vector a. Consider a small sphere 
around the point ( x,y,z ); then div a measures the flux of a 
leaving the point and crossing this sphere, in the limit, as the 
volume of the sphere shrinks to zero. We will show how the 
first term arises. 

Consider a small box, sides Sx, by, 8z, with the point (x,y,z) 
as one corner. Then the flux of a entering across the face perpen¬ 
dicular to the x-axis at the point (x,y,z) is a x dy.dz. The flux 
leaving the parallel face at x+8x is then 

( da x \ 

la x + — dxjdy.8z. 


Thus the net flux leaving the box is ( da x /dx).5x.5y.5z in the 
x-direction. This immediately gives the first term of div a as the 
volume SxdySz shrinks to zero, per unit volume. 

Finally, a vector is defined by 


, „ /. 8 . 8 , 8 \ 
curl a = v x a = i—+j — + k — 

\ 8x 8y 8zJ 

i j k 

= _i A jL 

8x’ 8y’ 8z 


:(ia x +ja y +ka z ) 


. / 8a z da \ ., 

-Ur 4 +)| 


< da x da z \ ,(8a y da x \ 
k 8z 8x J \dx 8y)' 


This vector is a measure of the work done in passing round the 
point ( x,y,z ) in a small loop perpendicular to the x-, y-, and z- 
axes in turn. 
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Reference should be made to books on vector methods for a 
full treatment. (D. E. Rutherford, Vector Methods, Oliver & 
Boyd.) 

It is frequently the case that we can define a potential function 
4>(x,y,z ) such that the force experienced by a test object can be 
written as a = —grad$. (For example, the electrostatic force on 
a test charge a distance x from a charge e is given by ejx 2 = 
(— djdx)(e/x), and ejx is the potential function in one dimension.) 
Gauss’s theorem shows that the flux of this force leaving a 
surface around a point is zero unless the surface encloses a 
“source”, for example, an electric charge, and results in the 
general law expressed by Poisson’s equation: 

— V . a = V 2 (j> = 4 np, 

where p is the source density. This gives Laplace’s equation: 

V 2 <£ = 0 

in a region without any sources. 

These two equations are of the utmost importance in engineer¬ 
ing, and describe the behaviour of such things as fluid flow, 
electrostatics, elasticity, and an immense variety of other problems. 

For reference, the form of V and of V 2 is given here in Cartesian 
coordinates, cylindrical polar coordinates, and spherical co¬ 
ordinates. 



(x, y, z) 
(r, 9, \j/) 
{r,9,z) 
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Note the different radius derivatives in V 2 for ( r,G,4'/) and 
(r,0,z) coordinates. 


SEPARATION OF VARIABLES 

Consider a two-dimensional form of the wave equation 
5 2 4> 1 d<j> 1 d 2 <l> 1 a 2 (j) 

dr 2 + r dr + r 2 SO 2 c 2 dt 2 


This could result from a problem with cylindrical symmetry, and 
with no variations in the z-direction. Suppose we may write 


<£ = *000(0)7X0. 

1 dR\ nrn /l d J 0\ 

f “r dr) + \ 7 2 d^r 


RQ d 2 T 
c 2 dt 2 ' 

„ , lfd 2 R Id R\ 1/1 d 2 0\ 1 d 2 T 

Therefore -^_ + - + _j , _ _. 

The l.h.s. is independent of t, so we may equate it to a constant, 

: 2 T = 0 (1) 


Therefore 


d 2 T 

dr r+ 


r 2 d 2 R r d R 
R dr 2 + R dr +M 


1 d 2 0 
= 0 'dO 2 ' 


Again, the l.h.s. is independent of 0, so we write 


( 2 ) 

( 3 ) 


The solution of (1) is T = Acoscnt, 

(3) is 0 = Bcoskd, 

(2) is a Bessel function of order k, 

R = CJ k (tir) + DY k (nr). 

Therefore <£ = (A cos cnt)(B cos k0 ) [CJ„ (nr) + DY k (nr)]. 
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It is frequently useful to separate the variables in a differential 
equation in this way. In the above example, suppose we were 
looking for solutions, periodic in time, in a cylinder of radius a, 
where $ is to be finite, single valued, and zero at the boundary 
at radius a. 

Then = A cos net cos kdJ k (nr), 

where n is to be chosen to give J k (na) = 0; this has first, second, 
..., mth zeros corresponding to values of n which can be found 
from tables, and we may describe the solution completely by the 
two integers k and n = m. 

The reader should not shy away from the perhaps unfamiliar 
words “Bessel functions”. Bessel functions are a set of wave-like 
functions which tend to appear in engineering and physics when 
a problem has cylindrical symmetry. They are tabulated (like sine 
and tangent) in books of tables. 

It is often the case that, as with a Fourier series, we can build 
up a solution to fit certain boundary conditions from a linear 
combination, that is, a sum of multiples, of such functions. These 
functions are called the eigenfunctions or modes, and the numbers 
corresponding to m or n, and k in the above are called eigen¬ 
values, or characteristic values of the problem. A simple example 
is a vibrating string; Pythagoras showed that the musical notes 
have a relationship described by the first, second, third,..., modes 
of vibration of a stretched string. Here the modes or eigen¬ 
functions take the form /Isinfax) and the eigenvalues m the value 
of n. 


FINITE DIFFERENCES 


The definition of the derivative f'(x) of the function f(x) is 


u f(x + 5x)-f(x) 


d/(x) 

dx 


In analogue work it is frequently the case that the differential 
equations describing the problem involve two or more independ¬ 
ent variables. Consider, for example, the flow of heat through a 
plate heated at one edge. The flow of heat is proportional to the 
temperature gradient, so this gradient must be written in an 
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approximate form by going back to the definition written above. 
The plate is divided up into a mesh of points, and the temperature 
gradient at a point written in terms of the difference of tempera¬ 
ture between two points and the distance between them. 

Thus, if x 0 , Xj, x 2 are three successive values of x, a distance h 
apart, 


d/i _ (x x )-/(xo) 
dx h 

_ f(x 2 )-f(x t) 
h 

/(x 2 )-/(x 0 ) 
2 h 


backward difference, 
forward difference, 
central difference. 


Hence 


_ 1 ( f(x 2 )- (x t ) /(x,)-/(x 0 ) \ 
dx 2 h\ h h ) 

_f(x 2 )+f(x 0 )-2/(x x ) 
fi 2 


It is readily seen that, in two dimensions, if f u f 2 , / 3 ,/ 4 are the 
values of /(x,^) at the points 1, 2, 3, 4 shown in Fig. 3, then 

= p (/i +fi +fi +A ~4/o) =/o 
And, in three dimensions, 

/S = p (/i +/ 2 +/3 +/4 +/s +fe - 6/o). 



Fig. 3. 
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A more sophisticated approach is to expand the function as a 
Taylor series about the point. Thus, in the lattice shown in Fig. 4 
we have 


- f _r n d fo“l 3 2 /o 

- f ° + “ ! Tx + v. IP' 


a\ d 3 f 0 
31 ~fa 3 + ' 


f 1 =f 0 -a 1 8 A + a l^_ a _l d Jfo + 

Jo 1 dx21 ox 2 3! 8x 3 + '" 


Therefore - a ^ + a 2 fi-{a l +a 2 )fo (a|-ai) d 3 f 

dx 2 \. a l a 2 (a l +a 2 ) 3(a 1 +a 2 )dx 3 ' 



Similarly, we may derive the ^-derivatives, and hence the value 
of V 2 / in finite difference form. The third derivative, which is 
neglected, will give an estimate of the error involved in the 
solution. 

This method is quickly extended to three dimensions or to 
other systems of coordinates. 


LAGRANGE COEFFICIENTS 

Suppose we are given the value of a function f(x) at the points 
x 0 , x t , x 2 ; then we may write 


(Y) (*-*iX*-* 2 ) 

(•*0 — *i)Oo — x 2 ) 


/(*o) + 


(x-x 0 )(x-x 2 ) 

( X 1 — X 0 )(Xj —x 2 ) 


/(*l) 


l (x _ xo)(x _ xl) 
(x 2 -x 0 )(x 2 -x 1 ) 


■f(x 2 ). 
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This has fitted a parabola through the three points. If more 
points had been taken, we could fit an n— 1 th order curve through 
the n points with the polynomial. 


/(*) = 


i 


(x-x 1 )(x-x 2 )...(x-x„) 

(x i -x 1 )(x i -x 2 )...(x i -x„) 


/(*.), 


where the term x—x t is omitted from the numerator, and 
(Xj-x ;) is omitted from the denominator of each ith term. 
That is, the coefficient of /(x ; ) is arranged to equal 1 if x = x i( 
and zero for any of the other points. 

By writing out the polynomial and differentiating, and then 
substituting the values of x, finite difference forms for the deriva¬ 
tive in terms of the points taken arise. The coefficients of each 
term are called Lagrange coefficients. A short table is given. 

Three points: 


2¥o= 3/0+4,-/o 

2hfl = - 0 + f 2 >+the order of 3rd derivative terms. 

2/1/2= /o— 4 /i + 3 /J 
and h 2 f" = / 0 - 2 /j + f 2 at all points. 


Four points: 

6h/o = — 1 l/o +18/1 — 9/ 2 + 2/jj 

6/1/1'=— 2/ 0 — 3/,+ 6/ 2 — / 3 I + order of 4th deriv- 

6 / 1/2 = fo~ 6 / 1 + 3 / 2 + 2/3 [ ative terms. 

6/1/3 = - 2/o+ 9/, —18/ 2 + II/3J 

/i 2 /o'= 2/ 0 - 5/,+4/ 2 - / 3 

/i 2 /r = /o- 2/, + f 2 

h 2 fi = /1- 2/2+ / 3 

/i 2 /3 = - /o+ 4/,- 5/ 2 + 2/ 3 


MATRICES 


Consider the simultaneous linear equations 

a u Xi + 012X2 + 0,3X3+ .. . + a ln x„ = c, 
o 2 i x i +a 2 2 x 2 + - ■ ■ +«2 „x n = c 2 


a„i X, +0„ 2 X 2 + . . . + a nn X n = C n 
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It is convenient to consider the quantity, called a matrix, written as 


A = (a y ) = 


The above equation is now written in terms of the matrix A and 
the vectors 


X = 

x l 

and c = 

c i 


x 2 


c 2 


x„ 


c„ 


or £ a t jXj = c h i = 1,2.n. 


That is, 




A can therefore be thought of as an operator relating x and c, or 
as an array of numbers, or as a kind of numerical name for the 
set of equations. 

If each element is replaced by the element mirroring it in the 
main diagonal, that is, a tJ is replaced by a jU the resulting matrix 
is called the transpose A' of the original A. 

The above defines a simple compact notation for the algebra 
necessary for dealing with simultaneous linear equations. It can 
be shown that the sum of two matrices is the matrix composed 
of the sum of the elements of the two matrices, thus 


(a i} )Hb iJ ) = {a ij +b i] ). 

If A = (aij) and B = {b i} ) are two matrices, the product C = B. A 
is derived as follows: 


The kth member of A. x is £ a kj x j- 
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So the /th member of B(Ax) is £ yb ik £ a kJ Xj 
That is, the product B. A = C = (c tJ ) is given by 

(cy) = ( !>*«*;) ■ 

Observe that BA # AB. 

The identity matrix is the matrix 

1 = 1 , 0,0 ... 0 
0 , 1,0 ... 0 

0 , 0,0 ... 0,1 

with l’s on the principle diagonal, and zeros everywhere else. 


DETERMINANTS 

Determinants are numbers which are connected with matrices. 
Consider the numbers defined by 



and so on, where the elements of the first row of the array are 
taken with alternating sign, multiplying a similar array, derived 
from the first by leaving out the first row, and the column in 
which the element occurs. 

These numbers are called determinants. To expand the deter¬ 
minant detA = det(a ;j ) = ||zz fj -|| the following method is used. 
For the terms in the expansion by the /th row, a u is multiplied 
by (- l) ,+J and then by the determinant obtained from the first 
by omitting the /th row and yth column. The terms obtained by 
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setting j= are then summed to give the value of the 

determinant. 

Suppose we have n simultaneous equations defined by Ax = c. 
Then, if an inverse A -1 of A exists such that A _1 .A = I, and 
it can be found, we may write 

A~‘Ax = A -1 c 
or, since A -1 . A = I, the identity matrix, 
x = A “ 1 . c. 

It can be shown that such an inverse matrix exists if and only if 
the determinant of A is non-zero. If det A = 0, the matrix A is 
said to be singular, and the solution has to be found by different 
means, depending on the nature of A. 

If we want to solve the equation Ax = 0, clearly x = 0 is a 
trivial solution. The condition that a non-zero solution exists is 
that detA = 0. 

EIGENVALUES AND EIGENFUNCTIONS 

Consider the relation A. x = y connecting the two vectors x and 
y. If any vectors exist such that A. x = Ax, where A is a scalar we 
may write 

(A — AI)x = 0. 

For non-trivial solutions, det(A-AI) = 0. 

This is a polynomial of degree «, which may be solved for its n 
roots. These roots A, are called the eigenvalues or the characteristic 
roots of the matrix A, and the corresponding vectors x t obtained 
by solving (A—AI)x f = 0 are called the eigenvectors. They have 
the property that any vector x can be written as a linear combina¬ 
tion of the eigenvectors, x = a 1 x 1 + a 2 x 2 + ...+a„x n , and that no 
eigenvector can be written as a linear sum of the other eigen¬ 
vectors. This is summed up mathematically by saying the 
eigenvectors span the space of vectors considered, and that they 
are orthogonal. 

A mention was made of eigenfunctions earlier. These two 
properties of spanning and orthogonality apply to any set of eigen¬ 
functions, in a modified form. The best-known set of eigen¬ 
functions is that used in expanding in Fourier series, sin«0 and 
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cos nd, where n goes from 1 to infinity. It can be proved that this 
set of functions span the set of all well-behaved functions with 
period 2n, and they satisfy the orthogonality relations 

J sin «0sin mOdO = 0 if n # m 

= n if n = m. 
f sin nO cos md dd = 0 


= n if n = m. 

These relations enable us to expand a periodic function f{6) in 
the form » ■» 

iA 0 +£ A n cos n9 +£ B n sin nQ, 


where A n 




f(9)cosn0dQ, 





/(0)sinn0d0. 
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chapter 2 


ELECTRONIC 

ANALOGUE 

EQUIPMENT 


THE AMPLIFIER 

THE D.C. AMPLIFIER 

The d.c. amplifier is the basic building block of the electronic 
analogue computer. It can be made to sum, multiply by con¬ 
stants, and integrate voltages, and to drive such things as multi¬ 
pliers, voltmeters and recorders in associated equipment. Much 
design effort has gone into the perfection of low noise, stable, 
drift-free amplifiers with band widths of up to 100 kc/s. Such an 
amplifier is a precision piece of equipment, very reliable, and 
useful in many fields other than analogue computation. 

Consider Fig. 5a; the round-backed triangle is the conventional 
symbol for a d.c. amplifier. We assume for the moment that it 
draws no current at its input A, which will in practice be the grid 
of its first valve, that it has an effectively infinite negative gain, 
and that it does not drift. If a small voltage V g exists at A, a “very 
large” negative voltage will be produced at the output B, and a 
current will flow along AR 0 B until V g becomes sensibly zero. This 
input point A is, therefore, called the virtual earth; for a practical 
computing amplifier with a gain of perhaps one million, it will 
never be more than about 01 mV from earth, which will, there¬ 
fore, give an output of 100V. 

The amplifier in Fig. 5a is fed by a voltage V, through a 
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resistor R u and has a feedback resistor R 0 from the output volt¬ 
age V 0 . Writing down a current balance for the virtual earth A, 
we have from Kirchhoff’s law. 


Therefore 


-v 0 = p Vl . 


Clearly, if voltages V u V 2 , V 3 ,..., feed the virtual earth through 
input resistors R u R 2 , R 3 ,..., we have (Fig. 5b) 


-HWMW 


or 

or 



-Vo = ^y t v b 


where y t = R 0 /Ri with the dimension of a pure number, is called 
the gain of the ith input, and the summations are taken over all 
inputs. 

An amplifier connected with resistors in this manner is therefore 
seen to act as a summer of the input voltages, with a sign in¬ 
version on the sum. The conventional symbol for a summing 
amplifier with gains y u y 2 , y 3 is shown in Fig. 5b. 

If we consider the general case of an amplifier with input 
impedances Z t and a feedback impedance of Z 0 , the current 
balance still holds, exactly as in the case of the summing ampli¬ 
fier. It can be seen that in this case we obtain 


where the summation is over all inputs. 

If the input impedances are resistances R h and the feedback 
impedance is a pure capacity C, this equation becomes that of an 
integrator, thus: 
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or *o = EVi —, 

i P 

where y ; = —. 

'* -R; C 

Since the transfer function between V 0 and the ith input is 

1 _y, 

pCRi p ’ 


which is that of an integrator, such an amplifier gives the sum of 
the integrals of the input voltages with respect to time. The 
constants , with dimensions of number divided by time, are 
called the gains of the integrator. 

Writing down the differential equations for one input of an 
integrating amplifier, we have, alternatively (Fig. 5c) 


—i +C- 


Therefore — V„ = —— f V, d t, 

Rx cj (0 1 

where t 0 and t l are the starting time and the time of observation. 

The conventional symbol for an integrating amplifier with 
integrator gains y 2 , y 3 , is shown in Fig. 5c. 

By putting a capacitor in the input and a resistor as the feed¬ 
back component of an amplifier, a transfer function of the form 
pCR is obtained. This corresponds to a pure differentiator. 
However, a differentiator accentuates any noise that may be 
present in the circuit, since its gain increases with frequency. 
Crudely, the differential of a spikey waveform is a much more 
spikey waveform. There is always some noise present at 50 c/s 
for a reason to be explained below, and the gain of a differentiator 
at 50 c/s is 100 times that at 0-5 c/s, which is a high frequency for 
an analogue problem, except in special applications. 

For this reason, pure differentiators are almost never successful 
in analogue circuits. All equations are transformed into terms 
of integration as the operation for solution of the differential 
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equations for the system variables. Occasionally a differentiation 
is unavoidable, and it must be approximated; this will be dealt 
with in a later chapter. 

BASIC RELAY CONTROLS OF THE AMPLIFIER 

We have shown how to use the amplifier for summing or 
integrating. In order to control the amplifier operation, relays 
are used to switch inputs to the amplifier. Figure 6 shows a 
typical arrangement. The common point of the input resistors. 



Fig. 6. 

called the summing junction, is connected to the virtual earth via 
two relay contacts. For the amplifier as shown, the problem 
clamp relay connects the virtual earth direct to earth. The output 
of the amplifier, therefore, remains constant, or clamped, at the 
voltage across its integrating feedback capacitor. If the relay is 
changed over to make contact to the summing junction, integra¬ 
tion starts. If the initial condition relay is changed over also, the 
amplifier gives as output a constant voltage determined by V it the 
initial conditions input. Releasing this relay then gives integra¬ 
tion from this voltage as an initial condition. 

Two further contacts enable the amplifier to be switched to give 
a capacitor or a resistor as feedback component, for its use as 
an integrator or as a summer. The operation of these contacts is 
self-explanatory. 

By means of more complex relay circuitry, computers are 
arranged to have various conditions or modes of operation. In 
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the “operate” or “compute” condition, summing junction and 
virtual earths are connected for all amplifiers. In the “hold”, 
“clamp” or “freeze” conditions, summers are in the operate con¬ 
dition, but the virtual earths of all integrators are earthed. In the 
“gain check” or “rate test” condition, all amplifiers are switched 
to the summing conditions, and the gains of integrators as well 
as summers can then be adjusted. “Set initial conditions” 
switches the integrators to charge their feedback capacitors to 
the required initial voltage, using an additional initial condition 
input. Finally, “standby” or “zero set” switches all amplifiers to 
the summing conditions, with only one input to each; those 
inputs are earthed, and arranged to give a high gain of perhaps 
1000. The amplifiers can then be trimmed to ensure that they 
give zero volts output for zero volts input. 

Exercises 

1. Modify the circuit of Fig. 6 to give a “zero set” facility. 

2. The facilities of “zero set”, “operate”, “gain check”, “hold” and 
“set initial conditions” are required for a computer with ten amplifiers. 
By operating switches, five out of the ten can be switched to be either 
summers or integrators, whilst the remainder are always summers. 
Devise relay circuits to give the required facilities, which are to be 
initiated by the operator closing a switch for the required condition. 


ERRORS OF THE SIMPLE AMPLIFIER 

The primary cause of error in d.c. computing amplifiers is drift. 
The circuit of a simple form of amplifier is given in Fig. 7. Suppose 
the input voltage V x is made zero and the potentiometer R v is 



Fig. 7. 
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adjusted until the output V 0 is made zero. If the amplifier is now 
left, the output will drift from zero because of such things as slow 
changes in the properties (the g m ) of the valves, alterations of the 
H.T. supply, or mains voltage drift resulting in heater current 
changes, temperature changes and so forth. The amplifier is 
particularly sensitive to changes affecting the first valve, a small 
alteration here affecting all subsequent voltage levels. Good 
design can reduce the drift of such a straight d.c. amplifier to 
values of 100/fiV per day; the techniques used are ready aged 
valves, balanced stages of “long-tailed pair” type, very good 
voltage regulation and temperature control. But the invention 
and perfection of the chopper stabilization method by E. A. 
Goldberg in the 1950’s vastly improved the performance of the 
d.c. amplifier, and made it into a reliable and accurate piece of 
equipment. 


THE CHOPPER STABILIZED AMPLIFIER 

An a.c. amplifier is basically free from drift, since it does not 
pass d.c. However, it can be used to give a drift-free d.c. amplifier 
in the following way. By using a chopper relay, the input of an 
a.c. amplifier is fed alternately from an input voltage V lf and 
from earth, at a rate of, say, 50 c/s. This a.c.signal is amplified, 
and then rectified and smoothed using a simple RC filter, and a 
second contact of the chopper relay to act as the rectifying 
element. 

The frequency response of this amplifier is limited by the 
chopper frequency and the smoothing filter, so that it is not, by 
itself, a practical amplifier. It must, therefore, be combined with 
a d.c. amplifier to increase its band width. 

Consider the combination shown in Fig. 8. The virtual earth 
of a d.c. amplifier is fed to a chopper relay, which feeds an a.c. 
amplifier with V g and zero alternately. The waveforms at A, B 
and C illustrate the actions of the a.c. amplifier and the chopper 
relay, to give voltages V A , V B and V c . 

The d.c. amplifier has a balanced differential input, with, in 
fact, two separate inputs. The output is given by the gain times 
the difference between the two input voltages. The d.c. amplifier 
is connected so that the smoothed output of the a.c. amplifier V c 
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now feeds its differential input stage, i.e. the amplifier now 
amplifies (V g —V c ) by —M, giving output V c = — M(V g —F c ). 
Now V c = — NV g , so the final output is 

V c = -M(N+l)V g . 

This gives a much greater gain, but it does more than that, for 
the drift is present in the d.c. stage only. If the total drift voltage 
of the d.c. stage is V d , referred back to the input as if it occurred 
there, we have 


V 0 = -M{V g +V d -(-NV g )}. 

Therefore V 0 = -M(N+ l)V g -MV d . 

Thus the gain for the drift term is reduced by a factor N +1 
compared to the total gain from the combined a.c. and d.c. 
amplifying stage. 

In practice typical values are N = 2000 and M - 30,000, giving 
a total d.c. gain 6 x 10 7 , the ratio of drift to signal being improved 
by a factor of 2000 over that of the d.c. stage alone. The circuit 
of a modern d.c. amplifier is given in Fig. 9. The input stage of 
the d.c. part is a long-tailed pair, the virtual earth being taken 
direct to a valve taking a very low grid current, and the smoothed 
a.c. amplifier output going to the grid of the other valve of the 
pair. The output stage is a cathode follower, which can give up 
to 10 mA without distortion in the range +100 to - 100 V at an 
output impedance of 1000 £1 or so. This impedance is reduced 
to the order of a milliohm by the high negative feedback intro¬ 
duced by the external feedback component. 

The presence of the synchronous chopper relay circuit intro¬ 
duces noise at the chopper frequency. This is normally the mains 
frequency, the chopper being driven from the heater supply, 
although 400 c/s choppers are used occasionally. 


ERROR ANALYSIS 

Consider an amplifier with voltage V g at its input grid and gain 
— M, drawing no grid current. 
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Then the current balance gives 


i 


V i-V g 

Zi 


= 0 , 


where Z 0 is the feedback impedance. 


Also, 

Therefore 


V 0 



'Z 0 MkoZi' 


For an amplifier with five 1 Mil inputs, two 100 kft inputs, 
two 500 kQ inputs and a 1 Mfl feedback resistor, which will give 
gains of 1, 10 and 2, the quantity 


10 1 

= 30(MQ) _1 . 


For an amplifier gain of M= 10 7 times, the error term 


Y° yi 

M' L Z t 


in the above equation for V 0 /Z 0 gives an error current of 
3 x 10 _6 Fo fiA with the resistor values given. Now the output V, 
of a perfect amplifier would be given by 


z 0 yz’ 


f I5> z yi 

f M -Z °yZi‘ 


Therefore the actual voltage for a 1 Mfi feedback resistor is 
related to the ideal voltage V, by 

V,= F 0 (l + 3xlO- 6 ), 


which is an error of 3 x 10~? per cent of the output voltage. 

If there is drift in the amplifier, it can be referred back to the 
input, and thought of as due to a voltage V d at the grid. 

Then V 0 = —M(V g + V d ), which leads to 


,= iZ; Z 0 Mfi?oZ, d = o Z t ' 
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Using the resistor values given above, the error due to drift is 
thus 30U <; Volts. A value of 100/iV for V d is fairly realistic, giving 
an error of 3 x IO -3 V, or 3 x 10 -3 per cent of a 100V output. 

Exercise 

Show that the error due to a grid current of 4 x 10 3 /jA is negligible 
compared with the drift error, for both a summer and an integrator 
with feedback capacity of 1 //F. 

In addition to the above causes of error, which directly involve 
the amplifier, the wiring forming the interconnections between 
amplifiers will have a leakage resistance of the order of 10 5 Mfi, 
giving leakage currents of 10~ 3 pA. This means that the leakage 
of an integrator will be 10~ 5 V/sec, or 1 per cent in lOOOsec. 

Exercise 

The output impedance of a cathode follower output stage of a com¬ 
puting amplifier is 50011. If the amplifier gain is 10 7 , show that the 
output impedance with a 1 MO feedback is about 2-5 x 10~ 4 n. 


errors of summers and integrators 


The value of the gain of the d.c. stage of an amplifier falls off 
with increasing frequency; to a fair approximation, we may write 


M(p) = 


M 

1 +pTi 


N(p ) 


N 

1 +PT 2 ’ 


where T, =3-3x10 3 sec, T 2 = 20 sec. 
Thus the total gain is 


M = M(l+lV){l+(p7yi + jV)} 

(l + pTiXl+pT,) 


= M, 


a+pr 3 ) 

(l + p7\)(l +pT 2 )’ 


where T 3 =0 01sec. 

Referring to the formula for finite gain error, the error is the 
inverse of the gain, so we obtain the error curve shown in Fig. 10a 
using the asymptote approximation. 


Exercise 

Show that the integrator error curve is as shown in Fig. 10b for an 
integrator with input resistors as given in the text, and a feedback 
capacitor of 1 n F. 
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Error, % 



Thus, at 100 r/s, or 16 c/s, a summer will have less than01 percent 
error, and the integrator error will be less than 0T per cent from 
5 x 10~ 4 r/s to 10 4 r/s. 

The noise introduced by the chopper will be of the order of 10 mV 
for a summer, 2 mV for an integrator. 


MULTIPLIERS 


THE COEFFICIENT POTENTIOMETERS 

The simplest form of multiplication of a voltage is that given 
by tapping off a fraction of it using a potentiometer. If, in Fig. 
11, R 0 = R u and the potentiometer slider was moved from one 
end of its track to the other, the output V 0 from the amplifier 
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would vary from nothing to - V u where V 1 is the voltage above the 
potentiometer. Thus, if we wished V 0 to be given by V 0 = —yV t , 
by putting a reference voltage of 100 V as V u and reading 
V 0 , the slider could be moved until the desired fraction was found. 


R 0 



Fig. 11. 


Of course, the resistor R x to the amplifier virtual earth will load 
the potentiometer, and the loading varies with slider position, so 
that the value of V 0 is not a linear function of position. But this 
is of no importance, since we are not interested in the position of 
the slider, only in the behaviour of the unit as a whole. 

The conventional symbol for a potentiometer is shown in Fig. 
11. It is standard practice to arrange the gains of an amplifier to 
be exact numbers such as 10, 1 or 0-1. To set the gain of an 
amplifier to some other value, a coefficient potentiometer is used, 
and set up in the loaded condition either using a high input 
impedance voltmeter on the slider or by reading the amplifier 
output, with a standard voltage at input. 

THE SERVO-MULTIPLIER 

If the slider of a potentiometer is driven mechanically to a 
position proportional to a voltage V x and a voltage V y is put 
across the potentiometer, its slider will be at a voltage pro¬ 
portional to the product V x V y This may be done by a position 
servo, as shown in Fig. 12. The amplifier takes V x and compares 
it with V s from a position potentiometer which has the reference 
voltage, usually 100 V, across it. The output of the amplifier feeds 
a drive amplifier, which drives a motor whose shaft carries the 
potentiometer sliders round the potentiometer tracks. The 
motor will turn until it has reduced the amplifier input V x -V a to 
zero. The potentiometers are then positioned at a proportion 



40 


ANALOGUE COMPUTING METHODS 


VJ100 of their travel, and the sliders of the ganged potentiometers 
will, when suitably loaded, reproduce this proportion of the 
voltage across them. Thus we obtain the product in the form 
V x V y l 100 from an output amplifier whose input resistor matches the 
inputs to the multiplier drive amplifier. Conventional symbols 
are shown in the figure. 

In practice the drive amplifier will have a velocity feedback 
input as well, to improve the response, and various shaping net¬ 
works may be added. A good multiplier can be made to have a 
static accuracy of about 0-5 per cent, and an amplitude response 
which is 3dB down at 400 r/s. However, the phase response is 



-loo (Vy) 

(-VkMoo) 


Fig. 12. 

2-3° in error at 10 r/s, which gives the effective limit, since a 2° 
phase lag in the tracking of the multiplier will give a large error 
in the product. This may be seen by considering a 100V peak 
sine-wave drive, with +100V on the first multiplying slave 
potentiometer. When the input sine wave reaches 0 V, the slider 
will be at a proportion of sin2° x 0 03 of the travel, giving an 
output of 3 V, or 3 per cent error. 

The step response of a multiplier of this type gives typical rise 
times of the order of 20msec; following a step input to the drive, 
the sliders would have attained 90 per cent of their final movement 
in about this time. 
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The advantages of the servo-multiplier compared with the many 
other methods which have been devised for multiplication are 
reasonable cost, reliability and the ease with which it can multiply 
several variables by the same input quantity, merely by ganging 
several potentiometers to the motor shaft. They are mechanically 
simple, and easily maintained, and, with proper protection 
circuits to prevent continuous rotation with the consequent wear 
on the potentiometer tracks, are almost foolproof. Their dis¬ 
advantages are that they are slow, of limited accuracy, and they 
introduce a certain amount of noise. 

Figure 16a shows circuits for four quadrant multiplication, 
division and square rooting. 

Exercise 

In the divide and square-root circuits given in Fig. 16 what are the 
input voltages to the amplifiers from the slave potentiometers? Devise 
circuits to give cube roots and other non-integral powers. 

CROSSED-FIELD MULTIPLIERS 

Special cathode ray tubes having perpendicular electric and 
magnetic field are produced for multiplication. An electron gun 
produces a beam of electrons with about l-5kV as the accelerating 
field. This beam is deflected in the 7-direction by the input V x 
on the X plates; they therefore have an 7-component of velocity, 
which interacts with an axial magnetic field proportional to V r 
This is produced by a coil wound round the tube carrying a 
current proportional to V y . The interaction gives a deflection in 
the 7-direction. This deflection is then reduced to nothing by a 
voltage applied to the Y plates, and it is easily shown that this 
voltage is proportional to the product V x V r The detection of the 
deflection can be done by masking half the screen of a phosphor¬ 
escent cathode ray tube and amplifying the output of a photocell 
to give the T-voltage in such a sense as to tend to bring the beam 
of electrons just to the edge of the mask. Special electrodes within 
the tube may also be used. Generally a pair of plates, one above 
the 7-axis and one below, are fitted within the tube, and an 
electronic servo brings the beam on to the axis by balancing the 
currents in the two plates. This method has the advantage that 
it does not depend on the behaviour of the phosphor. 
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Such multipliers are fast, their accuracy depending on the design 
of the electrodes (for low secondary electron current from the 
detecting plates to the Y plates, and correct axial field), and the 
amplifiers driving the coil and the Y plates. A typical perform¬ 
ance is 1° phase shift at 2 kc/s, 1 per cent accuracy in that range. 
They need careful screening from stray fields, and careful setting 
up. They have been used in computers and response analysers 
to give the reference and quadrature components of a system out¬ 
put when driven by a sine-wave input. 


HALL EFFECT MULTIPLIERS 

If a rectangular semiconductor wafer in the X-Y plane has 
steady current /flowing in its X-direction, and a magnetic field B 
exists perpendicular to the plane of the wafer, a small voltage 
will appear across it in the T-direction. This voltage is called 
the Hall voltage, after its discoverer, and is proportional to the 
product IB. More exactly, 

RJBx 10" 8 


where I is in amps, B is gauss, V in volts, t is the wafer thickness 
in cm, and R h is the Hall coefficient in cm 3 /coulomb. Using 
transistor amplifiers, several multipliers have been built of very 
similar performance to the crossed field cathode ray tube multi¬ 
plier already described. 


QUARTER-SQUARE MULTIPLIERS 
The quarter-square multiplier depends on the relation 

*y = i{(*+j0 2 -(*-j0 2 } 

for its operation. 

If the squares of the sum and the difference of the inputs V x 
and V y can be formed, a subtraction will give the product. Passive 
circuits can be constructed to perform the squaring from diode 
circuits which will be described later. It is then a simple matter 
to obtain the product voltage. 
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These multipliers can be made fast and accurate. For small 
inputs, the error can be large unless the circuit is modified some¬ 
what, due to the flatness of the square curve near the origin. This 
is overcome by the following method in the Reeves multiplier. 
We have 

4xy = [3(x+ y )+(x + y) 2 ] - [3(x -y)+(x- y) 2 ] - 6 y. 

The two quantities in the square brackets are produced by diode 
networks and the final addition gives the product. The advantage 
is that the two curves of 3(x+y)+(x+y) 2 and 3(x-y) + (x-y) 2 
cut the origin at an angle, thus greatly improving the performance 
at small signal levels. 

Triodes can be made to perform the squaring by using the fact 
that the anode current is approximately proportional to the anode 
voltage squared. This and similar methods are often used in 
repetitive computers. 

TIME-DIVISION MULTIPLIERS 

If a repetitive pulse has its width made proportional to one 
variable and its height proportional to the other, its area is a 
measure of the product. If this is done fast enough, and modern 
techniques make pulse repetition rates of 20 kc/s possible, an 
extremely accurate, fast multiplier results. The pulse may be tri¬ 
angular, or the mark-space ratio of a square pulse can be varied; 
the output is then smoothed in an RC circuit. There are various 
subtleties of the electronics in the practical design, which are be¬ 
yond the scope of the present book and will not be described here. 

The electronics of these multipliers must be watched, as they 
readily loose their adjustment. But they are widely used and are 
capable of accuracies of the order of 0 05 per cent full scale from 
0 to 10 kc/s. 

If V x is made to vary the pulse duration, this pulse train can be 
sent to several modulators to give pulse heights proportional to 
V yi , V y2 ,so that several products by V x can be formed as in 
the servo-potentiometer multiplication, to give products V x V y , 
v x v } , 2 , etc. 
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DIODE FUNCTION GENERATORS 
AND LIMITERS 

The gain of an amplifier can be varied as a function of the 
input voltage by arranging for diodes to switch in resistors in the 
feedback or input network. Consider the curve in Fig. 13a. It 
can be simulated crudely by the two straight lines shown, the 
transition or “break” from one to the other being at 50 V. The 
amplifier in Fig. 13b can be arranged to give an output voltage 



to) 



(b) 
Fig. 13. 
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defined by the two straight lines. If the slider of the potentiometer 
P is at -50 V and V x is just below 50V, at 49V, the point A is 
below zero, and the amplifier gain is R 0 IR X . If V x is at 51V, A is 
above zero and the diode conducts. The gain them becomes 
R 0 IR 1 + R 0 jR 2 , if r is small compared with R 2 . 

Thus, to simulate our function 

Ro = 25 

R x 50 

Ro «o_75 
^ + 50' 

So, putting R 0 = 1 Mft, we find R x = 2MI2, R 2 = 1 gives 
the required approximation to the curve. 

Clearly, the number of diode networks of this type can be in¬ 
creased. About ten is usually sufficient to give a simulation of 0-25 
per cent accuracy for a wide range of functions. 


R A R 



Ro 

Fig. 14. 


On a typical function generator of this type, both signs of the 
input voltage and + 100V and — 100V will appear on specified 
sockets or sets of tags. Also, sets of potentiometers for adjusting 
the break points where each straight line segment takes over 
from its predecessor, and potentiometers to adjust the gains to 
give the right slopes to those segments will be fitted. 

The limiter is a special form of function generator. An ampli¬ 
fier with a diode arranged as in Fig. 14 will have a gain of R 0 IRi 
so long as the output remains below V L . If the output increases 
further, the point A will then rise above zero, and the diode will 
conduct. The output then remains at V L provided R is small 
compared to R 0 . 
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A variety of limiting circuits are shown in Fig. 15. A dead space 
circuit, a back-lash and a friction circuit are shown. An under¬ 
standing of their operation is left as an exercise to the reader. 

RESOLVERS 

An important class of function generators is that for the genera¬ 
tion of the trigonometrical functions. Sines and cosines can be 
produced from specially shaped potentiometers or from linear 
potentiometers in which the slider is positioned by a linkage, 
swash plate or by other means so that the output voltage follows 
the sine of the angle represented by the input. Today, the most 
common form is the specially wound potentiometer. A continuous 
circular track is fed by + V and - V at opposite ends of a dia¬ 
meter, and two sliders mounted at right-angles on the centre 
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(a) 




Fig. 16. 


spindle give outputs of Fsinfl and VcosO, where 6 is the shaft 
angle. It is usual to drive the shaft by a position servo as in a 
servo-multiplier, with two sine-cosine and two linear potentio¬ 
meter tracks. This gives a unit, which is sometimes permanently 
connected, for polar to cartesian or cartesian to polar coordinate 
transformations. 

The conventional symbol is shown in Fig. 16b. 

Exercise 

Devise circuits to transform from polar to cartesian coordinates and 
vice versa. 

NOISE GENERATORS 

The normal technique for noise generation is to use a thyratron 
as the primary noise source. After amplification and level stabiliz¬ 
ation, a noise spectrum flat from about 30 c/s to about 3 kc/s can 
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be produced. A region, say from 300 to 500 cs is selected by a 
band pass filter and shifted by a chopper relay technique, working 
at 400 cs, to give the sum and difference frequencies of 400 cs with 
the noise waveform. After filtering off the high frequencies, 
Gaussian noise to about 1 per cent from 0 to 35 cs, can be obtained. 
This is a wide enough band for most analogue uses. 


OTHER FUNCTION GENERATORS 

A function generator of 1 per cent accuracy can be made by 
placing a shaped mask in front of a cathode ray tube. The light 
spot is made to follow the profile of the mask by means of a 
photocell and an electronic position servo as described in the 
crossed fields multiplier, and the resulting 7-deflection voltage 
gives the output for any 7-deflection. 

A related form of function generator of about 0T per cent 
accuracy can be made from an 7-7 plotter. The normal X-Y 
plotter is a device which positions a pen in its X- and 7-coordinates 
by two position servos, one driving a bar to and fro above the 
plotting table, the other a pen carriage along the bar. The pen 
can be replaced by a pick-up head which can sense the field about 
a conducting ink-line drawn on the paper. A radio-frequency 
current is then passed along the ink line. Thus, if the bar is 
positioned according to the 7-input, the pick-up head will be 
driven to the line and a suitable output from the 7-servo will 
give the appropriate output. A two-variable function generator 
can be made using resistance paper with silver ink lines as equi- 
potentials, and a pick-up probe mounted on the pen-carriage of 
the plotter. The output corresponding to the X and 7 inputs is 
then interpolated between the voltages of the two nearest ink lines. 

Exercises 

1. Calculate the resistors for a diode function generator to give 
Vo = (7,/10) 2 for V\ in the range 0-100V. Start the first line segment 
from a point on the input voltage axis to give the curve near the origin. 
What is the accuracy of the simulation? 

2. Calculate the resistors for a reciprocal function generator. Start 
the first line, not at the origin, but at the largest input voltage, and work 
back. 
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3. Show how to use a servo-multiplier to give squares and higher 
powers; show how a power series can be produced by this means. 

4. The gain of an amplifier is 

M ^~ M TT7r 

where T = 20 sec. It is used with suitable resistors to give a d.c. gain 
of ■/. What is the gain function y(p)? Derive y(p) for the more accurate 
value of M(p) given in the text. 

5. Use servo-resolvers to give the range, and azimuth and elevation 
angles to a target from the origin, given x-, y-, and z-coordinates. 
Simulate a correction for the earth’s curvature. 
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chapter 3 


SCALING 

AND 

APPLICATIONS 


INTRODUCTION 

We have discussed in the last chapter the electronic equipment 
that has been produced to perform the operations of addition, 
multiplication, division, integration and the production of non¬ 
linear functions. But we did not discuss how this equipment was 
to be used in the solution of the problems for which it was 
designed. It can quickly be seen that the three amplifiers (two 
integrators and a summer) in Fig. 17 will solve the equation 

and produce a solution of the form 

x = Asin((ot+(j)) 

where co 2 = y. 

Suppose that amplifier 1 has output x; then amplifier 2 gives 
— x, and amplifier 3 gives x, so that the output of amplifier 1 
is —yx, which is equal to x, and so the equation is solved. But 
this ad hoc method must be put on an organized basis if the 
solution of more general equations is to be attempted. The 
process is best illustrated by an example. 

Consider the equation describing the damped oscillatory motion 
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Fig. 17. 


of a weight on the end of a light spring, for displacement f(t) 
of its end 

mx + x+k[x -/(?)] = 0. 

Turn this equation round so that the highest derivative is alone 
on the l.h.s.: 

-mx = x+k[_x-/(t)l 

If the voltage signals representing the three terms on the r.h.s. are 
added into an integrator, its output will be x, which can be 
integrated again to give — x, and the r.h.s. can then be produced. 
Thus we can quickly draw the unsealed block diagram in Fig. 18. 



Fig. 18. 


The next step is to find the values of the gains shown in the 
block diagram appropriate to the equation constants considered. 
This is the process of scaling. The physical quantity x in the 
equation must be related to a voltage V x in the simulation. 


SCALING 

There are many different methods of scaling; they all depend 
on two different ways of looking at the same basic situation, i.e. 
(1) V x = b x . x volts 
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In (1), 6* has the dimensions of volts per physical unit; in (2), the 
physical quantity x is normalized and then multiplied by the factor 
100 V. Both methods have their disadvantages in practice, the first 
causes muddle over the dimensions, and is cumbersome in use, the 
second gets mixed up on the page when writing out the equations 
and when drawing out the block diagrams. A method is set out 
here which is fairly simple to use for working out gains, and is 
easy to write. Of the methods used by the author, it is the most 
convenient, and voltages appearing in the computing machine are 
quickly appreciated as the physical quantity they represent. But 
it is recognized that a group of computer users will have their own 
conventions and routines. What matters is not the notation, but 
the efficient use of the machines available, and the speed with 
which an operator, skilled or unskilled, can prepare accurate 
equations for the machine. 


SUMMER AND INTEGRATOR SCALING 

Suppose the constants of the physical system of a mass and a 
spring give the equation 

x+2-0x+16|>-/(0] = 0, 

where x is the displacement in inches, and t the time in seconds. 
Suppose, further, that we are only interested in oscillations and 
changes of x where a is less than 5 in. Then, for a computer with 
amplifiers linear in the range ± 100 V, we can work with (20x) volts 
to represent x in., so that when x is at its “full scale” value of 
5 in., the voltage representing x, that is (20x) volts, is 100V. 

We have chosen a scale for length measurements, but the 
equation involves rates of change also. It is therefore necessary 
to estimate the greatest rates of change of interest, so that the 
voltages representing x and x are within ± 100 V. Now it can be 
seen that the equation behaves similarly to the undamped 
equation 

x+16[x-/(0] = 0, 

and this equation has a natural frequency of ^/16 = 4rad/sec, so 
let us consider values of ratio of change equivalent to oscillation 
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of less than 5rad/sec; that is, consider as a first approximation 
that 

A = v4sin(5t). 

Therefore x=5Acos(50 

jc = —5 x 5^4 sin (5t). 

We have agreed that A, the amplitude, is less than 5in., so we 
represented x in. by (20x) volts. From those three relations it can 
be seen that x is less than 25 in./sec, so that we can represent 
x in./sec by (4x) volts. Again, x is less than 125in./sec 2 , so that 
we can represent x in./sec 2 by (4/5x) volts. 

We now have voltage variables which are proportional to the 
original physical variables. Also they remain within the computer 
voltage range for the expected variations in the values of the 
physical variables. We require the voltage equation relating these 
voltage variables to be identical to the physical equation relating 
the physical variables. 

We therefore take the original physical equation and replace 
its variables in turn by the voltage variables (writing them in 
brackets). To ensure that the equation as a whole is unaltered, 
every time a voltage variable is introduced in a term, any constant 
it introduces—within its brackets—must be removed from that 
term by division as a gain written without brackets. 

This is illustrated, stage by stage, below. 

Take each term in turn from the original equation, written out as 
x+20x+16x= 16/(0- 

The first term is x; our voltage variable is (fx), so the equation 
becomes 

(fx)x£+2-0x+16x = 16/(0- 

Similarly, for x, our voltage variable is (4x), so the equation 
becomes 

(0 8x)x 1-25+ 2-0(4x)x£+16x= 16/(0- 

Finally, x and f(t), which are lengths, must be scaled to (20x) 
and 20/(0 volts, giving 

1 -25(0-8x) + 0-5(4x) + 16(20x) x - 2 \ = 16[20/(0] x ^ 
or (0 8x) + 0-4(4x) + 0-64(20x) = 0-64[20/(0]. 
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The process may be summarized as follows. Take the physical 
equation. Replace the physical variables by the voltage variables 
and multiply each term by an appropriate factor, such that the 
original equation remains. This gives the voltage equation for the 
voltage analogue of the physical system, and it gives the constants 
for that voltage equation. These constants are the gains of the 
amplifiers involved. 

Rewriting the equation, we have 

—(0-8x) = 0-4(4x)+0-64(20x)—0-64[20/(f)]. 

We could, therefore, sum the terms on the right with gains of 
0-4 and 0-64 to obtain (0- 8x) as the amplifier output. 

If (0- 8x) is fed to an integrator with a gain of 5, the output will 
clearly be — (4x), and, by integrating this, again with a gain of 5, 
the final variable, (20x), is obtained. The circuit is completed by 
the external forcing function [20/(0], an d sign changing ampli¬ 
fiers as appropriate. 

Exercise 

Draw out the circuit of the above, checking each sign. 

It can easily be seen that the terms on the r.h.s. of the equation 
can be summed directly into the integrator which gives (4x). The 
gains become, therefore, 5x0-4 = 2-0 and 5 x 0-64 = 3-2. 

Exercise 

Draw out this circuit. How many amplifiers are saved ? 

TIME SCALING 

Suppose that the above circuit is set up and the gains set 
accurately. If the driving function is a step of 40 V, which corres¬ 
ponds to a step of 2in. in the position of the end of the spring, 
the system will oscillate with a frequency of about 4 r/s, the 
oscillation dying away with a time constant of 1 sec. 

This is too fast for the average recorder, though many special 
recorders are made for frequencies of this order. So everything 
must be slowed down for the purposes of the simulation. One 
way to do this is to replace the 1 fiF capacitors by 10/iF. Another 
is to divide the integrator gains by 10. Both methods replace the 
old rates of change in the system by new rates of change one-tenth 
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as fast, and so machine time /' and problem time t are related by 
the scaling equation t' = (10t)j i.e. lOsec on the machine is 1 
sec worth of “real” time on the physical system. 

So the final gains of the time-scaled circuit are 0-2, 0-32 and 0-5 
in the integrators, and 1 in the summer. These correspond to the 
gains of 2 0, 3-2 and 5 in the integrators and 1 in the summer in 
the original analogue. The final circuit is shown in Fig. 19. 



f=iot -[-20f(t>] 

Fig. 19. 


Time scaling can cause a lot of confusion. It is almost always 
easiest to work out everything in terms of “real” time, and then 
think of the capacitor rates of charge being reduced (i.e. the gains 
reduced or capacitors increased, in all integrators) if things change 
too fast for the simulation to follow. Conversely, the gains are 
increased if the physical system changes too slowly for conven¬ 
ience and accuracy on the simulation. 

Record charts from a machine might, typically, have a chart 
speed of 3in./min. If time scaling by a factor of 10 slower were 
done, as in the example, 3 in. of the chart would become a record 
of 6 sec of the problem solution. It is often convenient to time 
scale by factors of 60, so that minutes of the physical problem 
speed up to become seconds of the analogue, or slow down to 
become hours of the analogue solution. The resulting charts can 
then be marked appropriately in seconds or hours, though pro¬ 
duced in minutes. 

MULTIPLIER SCALING 

The output from the slider of a servo-multiplier potentiometer, 
matched suitably, is a voltage V z , given by 
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where V x and V y are the inputs. The factor 100, which has the 
dimensions of volts, appears from the mechanism of the multi¬ 
plier, since a multiplier multiplies by a proportion—the drive 
voltage V x is normalized as the function V x /100 of the potentio¬ 
meter travel. 

For example, suppose the coefficient of viscous damping in 
the term in x in the equation considered above is given by 
H = /i 0 (1-0 027’), where T is the temperature, in the range 
0-100°C. Then we must generate, in place of 0-8x^x(4x), the 
term 


0-8 


- (4i)—0-02: 


.2 (4 x)(T) 
4 (100) 


where the factor 0-8 comes from the scaling of the rest of the 
equation. The scaling factors are removed as gains at each 
replacement of a variable by its scaled value at each stage, and 
the original equation of the physical system preserved, exactly as 
given before. 

Similarly, for a division, V x divided by V y gives an output volt¬ 
age of 



so the factor 100 must be removed as a gain when writing the 
equation for the voltages. 

To summarize the scaling process, the operations involved are 
listed below. 

(1) Decide maximum values for each variable of the physical 
equations, including rates of change. 

(2) List physical variables and maximum values, and write 
down voltage variables in brackets. Each variable is multi¬ 
plied by an appropriate constant to give a voltage variable 
which always lies in the simulator voltage range, for the 
values of the physical variable considered. 

(3) Take each variable in turn, and substitute its voltage equi¬ 
valent in the system equations term by term. Each term 
must be divided by the appropriate constant so that the 
system equations are preserved unaltered. 



SCALING AND APPLICATIONS 


57 


(4) When two physical variables are multiplied, divide by the 
simulator voltage range when introducing the second volt¬ 
age variable. This gives the multiplier output with dimen¬ 
sions of volts, not (volts) 2 . Now multiply by the numerical 
value of the voltage range so that the term is unaltered in 
the equation. 

(5) Conversely, for a division, multiply by the voltage range, 
and remove this factor as a gain. 

(6) If time scaling is required, after the gains have been found, 
reduce integrator gains by a common factor for a simula¬ 
tion to behave slower than the physical system, and 
conversely. 


A WORKED EXAMPLE 

The uranium fuel elements in a nuclear reactor are encased in 
alloy cans which are pressed on to the elements. It is found that 
humps can appear on the cans, and that the angular size 4> and 
the normalized height x of the humps can be expressed by the 
equations , . 

(i) 

-2‘h. ( 2 ) 

where a is a parameter describing the elastic properties of the 
alloy can. 

We first construct the unsealed diagram of the simulation 
circuit of Fig. 20a. Suppose that by some means we have gener¬ 
ated (j> x in amplifier 1. Take 4>x and divide by x to get in 
amplifier 2. Multiply by a, and match the multiplier potentio¬ 
meter for constant load in amplifier 3. This gives <j>, which can 
then be integrated in amplifier 1 to give <j>x, solving equation (2). 
Now we know tj>x, so multiply it by x and match in amplifier 4, 
and then add a and <j>x 2 into an integrator, amplifier 5, to solve 
equation (1) for x. Two amplifiers are required to drive the x 
and a multipliers. Note that in this circuit, no voltage for dxjdt 
or (d/dt)(^x) is generated. So these variables need no voltage 
equivalents. 






Fig. 20a. 



(IO' ; o)—^6>=0 <50x)—^7>=0 

7 1 = 0-01 time scale times 6 yi — 0-06 

72 = 0 02536 y 2 = 01522 

y 3 = 0-009375 y 3 = 0-05625 

Initial condition: Amplifier 1, 39-27 V or <j> = in, x = 1 
Amplifier 5, 50 V or x = 1 


Fig. 20b. 





60 


ANALOGUE COMPUTING METHODS 



We have: 


a: 2 units take (100}x) volts i.e. (50a) 

<p 3 radians take (100^) volts i.e. (50 <j>) 

<t>.x 4 units take (lOO^x) volts i.e. ( 250 a ) 

a 0 01 units take (10 4 a) volts 


Then we have, scaling the equations, 




(100) '50x 10 4 


100 (25<£x)(50x) 1-875 


50 'df 
d 


(25 (j>x) = - 


25.50' (100) 

(10 4 a)(50<£) 


( 100 ) 


^ (50x) = 0-02536 —- — 0-009375(10 4 a) (4) 


( 100 ) 
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Now equation (3) is solved by amplifiers 1 and 3, whose output 


(— 504>). (10 4 a) 

( 100 ) 


(100 x5O0a) volts. 


Therefore the gain of amplifier 1 is =0 01. 
Similarly, the output of amplifier 4 is: 


(25<t>x).(50x) 

( 100 ) 


= (12-5<£x 2 ) volts. 


The gains of amplifier 5 which solves equation (4) are therefore 
y 2 = 0 02536 and y 3 = 0 009375. 

The final scaled circuit is shown in Fig. 20b, and, after time 
scaling by a factor of 6, gave the results shown in the charts of 
Fig. 20c. 

The analogue solution of these equations took a little over a 
working day. Its accuracy was limited by the computer available 
at the time and the presence of a dividing circuit. The results 
show that the hump in the can will grow without limit until the 
can bursts, or get narrower and flatten itself out according to the 
value of a. 


A SMALL PROBLEM 

A typical problem involving a non-linear differential equation 
is the determination of the closing time of the main stop and 
throttle valves in the steam pipes of a steam turbine. This valve 
is held open against a strong spring under normal conditions by 
oil pressure in a cylinder. If the turbine is tripped for safety 
reasons, the oil is rapidly drained from the cylinder, and the 
spring forces the valve shut against the steam pressure. The 
general arrangement is shown in Fig. 21. 

Now if x is the travel of the piston and its area is A and a is 
the area of the drain port which is opened to trip the stop valve 
shut, the oil flow Q is given by 

Q = Ax 
= ka^JP, 
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where A: is a constant and P is the oil pressure in the cylinder 
The force on the piston is given by 

F = AP 


a 2 x 2 x 2 

= A.-rr-i = 39-16-5 lb wt. 

k 2 a 2 a 2 

As the stop valve closes, therefore, the accelerating force is 
given by 

—mx — deadweight of piston 
—spring force 
+ oil force 
+friction force 

+ steam pressure force on the valve. 


A typical set of constants is given below: 


Oil-piston diameter 
Oil-piston travel 
System mass 
Deadweight force 
Friction 
Steam force 
Spring force: closed 
open 

Spring rate 


11 in. outside, 3-75 in. inside 
1-625 in. 

8211b 

7671b assisting opening 
4401b resisting closure 
19501b assisting opening 
34221b assisting closure 
56141b assisting closure 
13501b/in. 
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It is left as an exercise to show that the equation of motion of 
the system is 

2-115*= -39-16^ + 1350x4-265. 
a 


Taking maximum values, 

x m = 2in. work in (50x) volts, 

x m = 50 in./sec work in (2x) volts, 

x 2 = 2500 (in./sec) 2 work in (0-04x 2 ) volts, 


this gives the scaled equation 


d(2x)-: 


135° 2 , 265 

+ 2-U5‘50’ (5 X)+ 2-115 
- (0-04x 2 )+25-54(50x)+2-506(100). 


If the simulation is now slowed down by a factor of 500, the 
gains involved and the rate of stop valve movement come into the 
range of a feasible simulation. The circuit is given in Fig. 22. 

It is then a simple matter to take runs from the simulator for 
various values of drain valve area. Figure 23 shows a typical run 



A 0-0512 

B 0-005 

C According to orifice size 

D 0-05 


Fig. 22. 
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giving valve travel and velocity for a 1-75 in. diameter drain 
orifice. From such rims, the orifice giving optimum closing time 
without damage to the steam-valve seat due to excessive closing 
velocity on impact can be chosen. A laborious hand calculation 



is otherwise needed to work out step by step the movement of the 
valve. This type of calculation becomes a simple routine job for 
an analogue computer, with the possibility of useful answers for 
the designer by return of post if computer time is available. 

ESTIMATION OF MAXIMUM VALUES 

A problem is usually scaled wrongly once or twice before its 
completion. The maximum values must be estimated from con¬ 
siderations of the physical system and from the solutions or 
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estimated solutions of the equations. With many systems of 
importance this is comparatively easy, since some steady state 
will be known. But for non-linear systems the cut and try 
method is often the only way, and far quicker than involved 
analysis of the equations. 


METHODS OF SOLUTION 

It is in general an easy matter to draw out a circuit diagram 
for the simulation of a set of differential equations on the lines 
indicated in the two examples given. However, a more explicit 
classification of the methods used is in order. 

The “direct method” of solution is used to solve for functions 
of the independent variable, simulated as time, from the equations. 
For example, the coordinates of a projectile fired with velocity V 
at an angle of elevation 6 are 

x = Fcos d.t, 
y = Fsin O.t—^g.t 2 . 

So x and y can be generated direct from integrators producing 
VcosO.t and Fsin O.t from a constant voltage, and an integrator 
to give a ramp \g.t 2 for the integrator for y. 

The “indirect method” of solution is that used in the two 
examples given earlier. It was first proposed by Lord Kelvin 
when solving a second-order differential equation with his 
brother’s ball-and-disc integrator. He fed the input of his machine 
with a trial guess to the solution and recorded the output. This 
output was fed to the input as a second guess, and the solution 
produced by an iterative process on his open-loop machine. He 
then conceived the notion of connecting the output to the input, 
that is, of closing the loop, and to his delight the solution was 
obtained immediately. This step is quite obvious to those who 
have experience of feedback systems, but can be puzzling when 
met for the first time. The indirect method, Kelvin’s method, is 
the normal line of attack on a circuit. 

Often a problem variable can be generated by the solution of a 
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subsidiary or implicit equation. The dividing circuit of Fig. 24g 
solves an implicit equation 




r = 0- 


This method is often simpler and more accurate than generating 
a variable with a diode function generator. An example is the 
generation of the exponential y = e“*. Since dy/dt = — >>(djc/dr), 
the reader may quickly draw a circuit to give y knowing dx/dt, 
and the initial conditions. 

An extension of the implicit function method is the use of 
amplifiers without feedback. If, for example, (x, y) coordinates 
are to be converted to polars (r,9) then r = xcosO+ysind, and 
r = {x 2 +y 2 ) i . If these terms are generated with a servo-resolver 
and with multipliers, and added with opposite signs into an 
amplifier without feedback, the output of this amplifier may be 
taken as 9 and used to drive the servo-resolver. When the two values 
of r agree, the resolver stops turning, and a solution for 9 results. 

This technique has been used, for example, to simulate the 
behaviour of a boiler which first heated up the water as it passed 
through the system, then it boiled at a constant temperature 
dependent on the rest of the system. No explicit equation for the 
point where boiling started could be found, so the multipliers 
which fixed that point were driven from an amplifier with very 
little feedback (a gain of about 500) by comparing two values of 
the steam flow from the system. These values were equal if the 
equations were satisfied, and a solution was formed on the system 
by reducing the difference to a very small value. 

When a problem is scaled, it may happen that the gains of some 
of the integrators come out to be far too large for simulation. For 
example, the time constant for conduction through a tube between 
water and hot gas in a heat exchanger might be many times that 
of the gas itself. Consider the equation for the heat balance of 
the gas in a section of such a heat exchanger of length 1 unit. Then 
dG 

m g c g — = Mc g (Gj-G 0 )- H{G- T ), 

where m g c g is the thermal capacity of the gas, M is the gas mass 
flow, H is a constant, G is the mean gas temperature between 
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inlet G; and outlet G 0 , and T is the tube temperature. If the gas 
time constant (its thermal capacity) is very small, that is the gains 
of an integrator giving G would be very large, we may write, 
approximately, 

Mc g (G i - G a ) = H(G - T). 

So the gas temperature can be generated with a summer. As 
another example, when considering the changes in power from 
a nuclear reactor due to the daily variation of load, the time 
constants of the fuel rod and can, which are about 2 sec, can be 
neglected. 


SOME STANDARD CIRCUITS 

Figure 24 shows some standard circuits which may be regarded 
as the blocks from which to form a simulation circuit. The 
summer, integrator and lag are straightforward. Figure 24d is a 
lag-lead circuit which can be used as an approximate differentiator 
with V 2 V 3 = 7iV 4- The use of positive feedback to achieve a 
perfect differentiator when y = 1 is shown in Fig. 24f. The divid¬ 
ing circuit shown gives a square root if the multiplier drive is the 
output. 

It is sometimes necessary to simulate a transport delay. These 
give an equal time lag to signals of all frequency, with no loss of 
amplitude, and have transfer function given by 

V 0 = e~ pT V h where T is the delay time. 

This cannot be simulated exactly with summers and integrators. 
Some units have been designed which produce a delay with 
magnetic tape by setting the read and write heads at a distance 
apart, or with curve following devices on the same principle. But 
the most common approximations are derived from the Pade 
expansion of e“ pT and are: 

__„ r 1—0-5 Tp delay = T sec, maximum phase 

~ 1 +0-5 Tp shift 45° (5 per cent error). 


_ pT ~ 12 — 6pT + p 2 T 2 delay T sec, maximum phase 
~ H + 6pT + p 2 T 2 shift 155° (5 per cent error). 






SCALING AND APPLICATIONS 


69 


This transfer function corresponds to the equation 

x(t) 12 pT 

y(t) 12 + 6 pT-pT 2 ’ 

12 pT 

bUt "12 + 6 pT-pT 2 ' y ~ V; 

therefore x = y + v; 

6i> I2v 12 

therefore -—.y. 

The circuits for these approximations are given in Fig. 24. 


Exercises 

1. Give the gains and circuit of an oscillator to produce an a.c. 
signal of variable amplitude with frequencies from 1 c/s down to 1 cycle 
in 10 min, using standard integrators and summers. How does the 
leakage resistance of the integrator capacitors affect the performance 
of the oscillator? How could this resistance be estimated? If possible, 
devise a circuit to compensate for this loss continuously, and test the 
circuits on a computer. 

2. Give a circuit which will give the in phase and quadrature terms 
of the output of a system driven by a perfect oscillator by generating 
the first two terms of the Fourier series expansion of the output. 

3. Give circuits to transform polar coordinates to cartesian co¬ 
ordinates and vice versa, using sine-cos resolvers or function genera¬ 
tors or other means. 

4. Investigate the equation x +2lco„x + co n 2 x — fit) on an ana¬ 
logue computer with sinusoidal and step inputs. Show the phenomena 
of resonance, critical damping, etc. 

5. Write down the equation of motion of a simple pendulum of 
length a. Draw scaled circuits assuming sin 6 = 0, sin 9 = 9 — J0 3 , 
and then solve the exact equation by transforming to (x, y ) coordin¬ 
ates and using i = a cos 6.6, y — —a sin 6.6. How does the 
period vary with amplitude? 

6. Show that the simultaneous linear equations 

’LaaXj - a = 0 

where i, j take the values 1,2may be solved by (1) generating 
Xi from the /th equation in a summer, (2) solving the equations 
TanXj — a = dxijdt with integrators. 





Fig. 25. 
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It may be shown that a solution is not necessarily obtained if the 
matrix of coefficients (an) is not positive definite. A solution is always 
obtained from method (2) with 




d y* 
d /’ 


7. By sketching circuits, find the roots of the equation 
P n (x) = x n + a n - ix”- 1 +... + aix + a 0 = 0 

by writing 

Pn(x) = {x[x(x . . . (x + an- 1 ) . . . + «3> + « 2 ] + «l} + « 0 . 
Generate P n (x) with multipliers and find the roots by varying the 
input. Hence scale and solve 

P 4 (x) = x* + x* - I9x^ + 11* + 30 = 0 
(roots -5, -1, 2, 3). 

It can be shown that, to obtain useful results when solving a poly¬ 
nomial, much higher precision is required in the working than is 
given in the resulting root (see Modern Computing Methods , National 
Physical Library, pp. 59-60). It is therefore normally impracticable 
to solve matrices or polynomials above the third or fourth order on 
analogue computers. 

8. Draw circuit diagrams for the equations: 

x - c(l - i* 2 )* + x = 0 (Rayleigh’s equation), 

X + ki + (A + Bx 2 )x = F(t) (Dufferin’s equation), 
x - s(l - x 2 )x + x = 0 (Van der Pol’s equation), 

fix + tx + (fi - n 2 )x = 0 (Bessel’s equation). 

9. Derive the circuits and transfer functions shown in Fig. 25. 

10. A servo-controller has transfer function 

do K (1 + 30p) (1 + 13/>) 

T t 1 p(\ +p) ‘(1 + l-3p)' 

Draw a simulation circuit, and scale for 

\T e \ < 25°C, 0 < 6o < 1. 

Add a limit so that the rate of change of output has a limit given by 

\p6o\ < 1-6. 
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COMPUTERS AND 

ASSOCIATED 

EQUIPMENT 


GENERAL DESIGN OF A COMPUTER 

A computer user, when studying a machine for the first time, 
will discuss the specification of the amplifiers and multipliers, and 
the provision made for other equipment such as function gener¬ 
ators, relay amplifiers, etc. A user wants a simple, clear patch¬ 
board, accurate means of gain setting, and simple control over 
the computer. Amplifier gains should be standard, or standard 
except for a few selected amplifiers which can have special com¬ 
ponents wired in. An adequate number of amplifiers must be 
available for the type of problem the user is solving; and problems 
have a habit of growing as additional complications are added. 
Maintenance and servicing of the machine should be simple and 
rapid. And, of course, the price must be low. 

PATCHBOARD 

A typical modern computer of more than about thirty amplifiers 
will have a central patching bay, with a removable patchboard for 
making the amplifier interconnections. Problems can be set up 
away from the machine on the patchboard, and, once tested, can 
be stored for future use. The patchboard usually groups a set of 
amplifiers together, with from five to ten inputs and four or five 
output sockets for each amplifier. These groups will be separated 
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by potentiometers, sets of commoned sockets joined together, 
reference voltage sockets, and multiplier potentiometers or other 
multiplier sockets. Other sockets will give lines to recorders, 
switches, relay contacts, diodes, or to outlet sockets which may 
be joined to another computer working in conjunction with the 
first. The sockets are usually marked and colour coded to 
simplify patching. 

CONTROL 

The control panel of a computer will give operating conditions 
for “compute”, “hold”, “stand-by”, “gain set” and so forth. 
Control is via push-buttons or switches, operating relay interlock 
systems which in turn work relays on the amplifiers to give the 
appropriate computer conditions. Lights indicate the condition 
the computer is in at any time. Some computers can be operated 
repetitively, that is, the problem set on the machine is solved once 
every second, or every minute, or every time a certain variable 
reaches a given level. Another repetitive mode of solution resets 
the machine at, say, fifty times a second, and the solution is 
displayed on an oscilloscope. Usually two or more machines can 
be coupled together and controlled from one central desk. In this 
way, a basic computer of say, fifty amplifiers, can be used with 
two others to give a computer of 150 amplifiers, 

MEASUREMENTS 

A high impedance digital voltmeter giving a reading of voltage 
to four figures with sign is usually incorporated in a computer. 
Alternatively a sensitive voltmeter can be fitted. By balancing for 
zero voltage difference between the slider of a calibrated potentio¬ 
meter or a decade resistance box, and the unknown voltage, 
accurate measurements can be made without loading that un¬ 
known voltage. For gain setting, either a known voltage can be 
connected to any input whilst all other inputs are earthed, and the 
output of the amplifier monitored as the gain is set to the required 
value with the potentiometer, or the slider of the potentiometer 
can be monitored, taking care not to load it with the voltmeter, 
and set to give the required fraction of the input voltage. In the 
first case, relays make all integrators into summers in the gain set 
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condition of the computer. Thus the gains of integrators are set 
up exactly as are summers. If the potentiometer sliders can be 
read direct, this is not necessary, as the voltage read at the loaded 
slider gives the gain directly. 

RECORDS AND READOUT 

The normal method of recording the results of a computer run 
is on strip chart recorders. A moving coil meter movement drives 
a pen or a hot wire across a moving paper chart to give a con¬ 
tinuous record of the variable. Alternatively, a pen is positioned 
over a sheet of graph paper in the X- and Y-directions by two 
servos driven by two machine variables, or arranged to give a 
record with the Y-direction driven by a time base. Some com¬ 
puters have high speed automatic print-out of the voltages on all 
amplifiers. This printed record can be kept with the chart records 
of the runs, or used as a guide for checking the accuracy of the 
machine and the gain setting. 

OVERLOAD INDICATIONS AND ALARMS 

In order to ensure that the computer is operating within its 
voltage range, a buzzer or flashing light is arranged to show when 
any amplifier is outside its operating range. A convenient method 
is to trigger a neon light from the output of the a.c. channel of the 
amplifier, which produces large pulses when the virtual earth is 
lost due to excessive input voltage or too high output current, or 
other misuse of the amplifier. Similar alarms show when a 
multiplier is being driven off the end of the potentiometer track 
or outside its operating range. The power supplies and voltage 
reference supplies will normally have some form of overload pro¬ 
tection, and some form of automatic sequencing will ensure that 
the computer is turned on and off in the right manner. 

MAINTENANCE 

The maintenance of a computer should not normally worry the 
operator. A routine maintenance scheme, and a full-time main¬ 
tenance section are necessary for large installations, and for a 
smaller computer, faulty equipment is often maintained by the 
manufacturers. 
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SOME TYPICAL COMPUTERS 

THE PACE SYSTEMS TR-10 

This computer is a compact transistorized desk instrument for 
the solution of small problems involving up to twenty amplifiers. 
The machine unit is 10 V, which leads to an accuracy of 0T per 
cent with the high quality resistors and capacitors used in the 
machine. Twenty coefficient potentiometers are supplied, and 
multiplication is performed with diode quarter-square networks. 
The amplifiers are patched into the circuit required using plug-in 
resistors and patchleads to capacitor units and potentiometers. 
The arrangement of the computer can be seen in Fig. 26; the 
amplifiers are at the foot of the vertical panel with the capacitor 
units, the diode multipliers and function generators, reference 
voltage points, comparitor relays, switches and common points 
appearing above them in plug-in units. At the top of the panel 
are the potentiometers and their patching sockets. The control 
panel is extremely simply arranged. A three-position switch gives 
reset, hold and operate conditions of the computer. Gain setting 
is done with an accurate graduated potentiometer and a multi¬ 
range meter arranged in a null circuit. A multi-position switch 
gives the output of each amplifier and the values of the power and 
reference voltage supplies. Any amplifier overload is indicated by 
a light on the panel. The computer is turned on by a single main 
switch, and uses about 60 W. 

THE SOLARTRON SC-30 (FIG. 27) 

This machine incorporates all the features of many much larger 
computers. It is a thirty-amplifier unit, which can be used in 
conjunction with other computers to form a much larger machine. 
Twenty of the amplifiers are sum/integrate units with five inputs, 
values 3x1 MQ, 2x01 MQ (0-1 per cent), and four feedback 
components of 1MO, 01MQ, 1/rF, and 0-1 fiF (0T per cent). 
The remaining ten amplifiers are summers with the same range 
of gains. Sixty coefficient potentiometers are provided. Four 
servo-multipliers or various combinations of mark-space and 
height units of an accurate electronic multiplier can be fitted. 











Fig. 27. 

computer to potentiometer set, problem check, compute and hold 
conditions, with indicating lights showing the computer state. The 
computer can be arranged to work in a repetitive mode, either 
resetting itself or resetting at fixed intervals. It can also be 
arranged to work with one part of the computer controlling the 
other, via the reset relay lines which appear on the patchboard, 
or to give a voltage scale change by a decade. Five amplifiers can 
be used in a non-standard manner, with any circuit wired round 
them on the special component trays provided. The whole is 
housed in a desk unit 80 in. x 36 in. high x 38 in. deep, with power 
supplies and adequate desk surface. 


The computer has a removable patchboard with amplifiers, multi¬ 
pliers and standard limiters, function generators, relay compari- 
tors, and 100 V reference supplies simply laid out and colour coded 
on the patchboard. A five-digit digital voltmeter is provided for 
gain setting and voltage measurement, but a meter and nulling 
potentiometer can be used if preferred. Push-buttons set the 
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THE PACE SYSTEMS 231-R (FIG. 28) 

Electronic associates have pioneered high accuracy analogue 
computers, and, from their standard range, a computer can be 
built up to fit the prospective users’ requirements. A typical Pace 
computer might have 100 amplifiers, 150 potentiometers, 20 or 30 
multipliers or resolver units, 20 diode function generators and the 
associated racks, power packs, digital voltmeter and display 
equipment. The computer has an all-metal removable patch¬ 
board, the interconnections being made by shielded patching 
wires which considerably reduce the circuit noise. The unit 
might have automatic voltage print-out, X- Y plotters, a mainten¬ 
ance bench and equipment, and perhaps two-variable function 
generators, electronic multipliers, and other equipment from the 
manufacturers range. 

The unit can be extended by the use of the interconnecting 
trunk wiring to give a much larger unit, controlled from the 
central desk. A wide range of amplifiers, multipliers, servo- 
multipliers and resolvers, electronic multipliers, function gener¬ 
ators, X- Y plotters and recorders is available. A notable feature 
is the automatic potentiometer setting system, ADIOS (Automatic 
Digital Input-Output System). A paper tape can be prepared, 
which will set all the gains of a problem in a matter of minutes by 
means of servo-motors on the potentiometers. With its aid, a 
problem can be taken off a machine and another put on, ready 
for production, in 20 min or so. This overcomes the major 
inflexibility of most analogue machines, since the inevitable 
machine idle time which occurs during a problem solution can be 
used for other problems, and the time booking system usual on 
digital machines can be used. 

The standard checking, gain set, computing and problem hold 
facilities are available, using a push-button control system, for 
real time solutions. Alternatively, the machine can be used 
repetitively at from 10 to 50 c/s of problem solution, using a 
separate display. The cabinet containing all computing elements 
is thermostatically controlled to reduce any possible drift resulting 
from temperature changes. The units slide out of this cabinet for 
maintenance, and replacement units can be plugged in at once. 
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The high standard of the Pace system is indicated by the 
manufacturers’ claim that they have supplied 70 per cent of the 
world’s analogue installations. 


THE ENGLISH ELECTRIC SATURN (FIG. 29) 

Saturn was designed for the simulation of nuclear power station 
plant with the simulation of a complete nuclear station in mind. 
It has six computers of 252 amplifiers which can be used 
together to form a single unit of 1512 amplifiers, and a large 
number of multipliers and non-linear devices. It is permanently 
housed in the end of a laboratory building, with the amplifier and 
component racks together with the power supplies and the main¬ 
tenance facilities upstairs, and the computer control room with 
the six individual control panels, recorders and patchboard units 
on the ground floor below. A single computer comprises 252 
amplifiers divided into twelve racks of twenty-one amplifiers, 
and 600 potentiometers, fifty per rack. High-performance 
servo-multipliers are used, there being three available per rack. 

A detachable patchboard, about 2ft square, is subdivided into 
twelve sections, one section of which is shown in Fig. 30. 

The functions of the various amplifiers, and their gains, are 
shown on page 83. Thirteen of the twenty-one amplifiers can 
be used as integrators or summers, three are for use with the 
multipliers, two are for use as function generators, comparison 
relays, or for other special uses, and the remaining three are 
summers. 

The multipliers and potentiometers are indicated on the dia¬ 
gram. The control panel has a subsidiary patchboard connected 
to the main patchboard for use with the switches and voltage- 
source potentiometers on the panel. The outputs of every amplifier 
are brought to a telephone-type jacking board between the ampli¬ 
fier and the patchboard socket of its output; this enables voltages 
to be injected using jacks at any point of the circuit when setting 
up or during runs on the machine. Control push-buttons on the 
desk give clamp, operate, gain set, operate in one-tenth time, set 
initial conditions, zero set, and other computer conditions. The 
gain-setting potentiometers are in a short corridor adjacent to 
the control desk, and gains are set using the computer digital 




Fig. 29. 
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Saturn patchboard layout for one rack of equipment. 




COMPUTERS AND ASSOCIATED EQUIPMENT 


Amplifier 

number 


01 

02 

03 

04 

05 

06 

07 


09 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 


Multiplier drive amplifier. 

Summer gains of 1, 10, 10. 

Summer gains of 1, 10, 10. 

Summer gains of 1, 10, 10. 

Summer/integrator gains of 1, 1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Function generator or special purpose amplifier. 

Multiplier drive amplifier. 

Summer/integrator gains of 1,1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1,1, 1,10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1,1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 5x1; 2x10; 2x2 (for 
multiplier) Initial Condition; or one-tenth of these 
values at Integrator. 

Summer/integrator gains of 5x1; 2x10; 2x2 (for 
multiplier) Initial Condition. 

Function generator or special purpose amplifier. 

Multiplier drive amplifier. 

Summer/integrator gains of 1, 1,1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1,1, 1, 10, 2 (for multiplier) 
Initial Condition. 

Summer/integrator gains of 1, 1, 1,10, 2 (for multiplier) 
Initial Condition. 

| — | ^ | +~[ : - 100 V, signal earth, +100 V sockets. 


A0 to A9, B0 to B9 ... E0 to E9 are the references of potentio¬ 
meters. 

Row D gives lines to the control panel auxiliary patchboard. 
Diodes for limiters are shown in rows C and D. 
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voltmeter. Covers prevent subsequent movement of the potentio¬ 
meter dials by mischance. A panel shows when any amplifier has 
overloaded, and another panel below indicates whether the multi¬ 
pliers are operating correctly. Twelve recorders are mounted on 
the panels above the desk, and other recorders can be used if 
desired. 

The whole computer is serviced by maintenance staff, who also 
turn on and check the computer every morning. 

Operating experience on this machine has shown that machine 
outage can be as low as 5 per cent. Typical usage of the machine 
is, however, 40 per cent production time, 40 per cent problem 
changing time, 10 per cent idle time, 3 per cent wasted due to 
incorrect usage, and 7 per cent fault time. 


THE ROYAL AIRCRAFT ESTABLISHMENT TRIDAC 

Tridac was designed for the simulation of aircraft and missile 
systems, although it can be used for any other type of analogue 
problem. It was among the first large computer systems designed 
in this country, and is one of the largest in the world. However, 
its design and performance seem in many respects out of date by 
modern standards, and it is chiefly of historical interest here. 

The computer has about 650 amplifiers, operating to a machine 
unit of 30 V, and a number of servo-multipliers and resolvers. It 
is housed in a special two-storey building at R.A.E., Farnborough. 
The amplifiers and other electronic units are mounted in cabinets 
which divide up into eleven “rafts” which vary in detail. Each 
raft has a patching unit, which is connected to the central control 
and patching desk. The multipliers and resolvers which are 
hydraulically operated are in an annexe with their power supplies. 

The whole computer can be controlled in subsections, and this 
control can be patched on the control desk. Various repetitive 
computing modes are available, and standard X-Y recorders, 
oscilloscopes and strip recorders give the output of the computer. 
Overload alarms and indications are given at the desk and on the 
rafts. A regular maintenance system ensures maximum reliability. 

The basic design of the computer started in about 1950, and it 
was being operated by about 1954. Since then it has been engaged 
on a wide variety of missile and aircraft problems. Modifications 
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to the original design of the computer have made it more versatile; 
it was originally built as permanently connected units for the 
simulation of the aerodynamics, the axis transformation, radio 
system and so forth. This type of layout was widely used in early 
computers, and has been almost entirely replaced by the more 
flexible arrangement of central patching with every amplifier 
available for any use. 

There are many other commercial computers, and there is not 
space to describe more than the few of them mentioned. The ones 
mentioned are fairly well known, typical, and up to date, but are 
not specifically recommended. 

Exercise 

As an exercise for the student, or as a project for a group, the design 
of a computer of, say fifty amplifiers is recommended. The aim should 
be to obtain maximum accuracy and the greatest ease and efficiency 
of operation. Assume that standard high-performance amplifiers and 
multipliers, and the required power supplies are available. 

The design will divide up into firstly, the layout of the operators 
controls, the patchboard, the number of multipliers and potentiometers, 
and the number and types of non-linear devices required, and, secondly, 
the relay circuits and switches, the cable runs and wiring, the choice of 
components and hardware to realize the requirements, and the layout 
of the parts of the computer to match the specifications decided in the 
first part. The internal design of the computer should aim at ease of 
maintenance, and simplicity and standardization of wiring. The large 
number of wiring joints must be made in such a manner that faults in 
wiring can be easily located and rectified. Such problems as lead 
resistance, possible earth loops and other sources of noise, temperature 
stability and reliability of components must be borne in mind. These 
and the many other problems of design give scope for a great deal of 
imagination and inventiveness in their solution. 
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chapter 5 


APPLICATION OF 
LARGE ANALOGUE 
COMPUTERS 


A NUCLEAR REACTOR SIMULATION— 

AN EXAMPLE 

In order to describe the approach required when a large number 
of amplifiers is needed to simulate a system, an example is worked 
through as if it were to go on to a computer. A nuclear reactor 
simulation has been chosen, since, although reactor simulation 
studies are but a small part of the total field of application of 
analogue computers, the method of preparation of data and 
equations and the procedure on the machine is typical. 


THE PHYSICAL SYSTEM 

In the nuclear power stations of the Calder Hall type, heat is 
generated in natural uranium rods by nuclear fission. The rods 
are about 1 in. in diameter and 3 ft long, and are enclosed in 
magnesium alloy cans. The cans have various fin arrangements 
to improve their heat transfer properties to the coolant carbon 
dioxide gas. These cans of uranium fuel are stacked, one above 
the other, in vertical channels of about 4 in. diameter cut through 
graphite blocks, which are about 8 in. square and 3 ft long. The 
graphite blocks are built into a cylindrical stack, about 30 ft high 
and 40 ft in diameter requiring about 40,000 blocks. The whole is 
86 
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enclosed in a spherical or cylindrical pressure vessel of thick steel 
sheet, which is connected via ducts to separate boiler vessels. The 
pressure vessel is filled with carbon dioxide at about 250 lb/in 2 
which is driven through the graphite channels and round to the 
boilers by axial or centrifugal gas blowers in the boiler ducts. 
The gas, which is heated by the fuel elements in the channels, is 
used to generate steam in the boilers, which in turn is used to 
generate electricity in conventional turbo-generator sets. 

The reactor generates heat by the process of fission. A uranium 
atom splits, emitting very fast moving neutrons. These fast 
neutrons, on average about 2 - 5 per fission, are very unlikely to 
cause another fission in natural uranium; they pass into the 
graphite moderator blocks and are slowed down to the thermal 
energy appropriate to the temperature of the graphite. These 
slow neutrons are very much more likely to produce fission, but 
only about 40 out of every 100 fast neutrons survive the slowing 
down process and cause fission. 

When the number of fissions occurring at any instant is such 
that an equal number of fissions is caused as a second generation 
by the neutrons resulting from the first fissions, the reactor is self- 
sustaining, and the average rate of fission and the average neutron 
population remains constant. The reactor is then said to be 
“critical”. If neutron-absorbing material is removed from the 
reactor core, the neutron population will increase; this increase 
will cause a further increase, and the rate at which fission occurs 
will increase more and more rapidly in an exponential manner. 
The reactor is then said to be “supercritical”. Conversely, if 
absorber is inserted, the number of neutrons present and the 
fission rate will fall, and the reactor is subcritical. The number of 
neutrons available for fission from the preceding generation 
measures the reactivity. The natural reactivity of the reactor will 
be about 101 neutrons per neutron, and this 1 per cent super¬ 
criticality is controlled by the absorbing control rods and other 
effects. 

The problem is to simulate the dynamic behaviour of the 
reactor, its control with the neutron absorbing control rods, and 
its cooling and safety as regards the gas flow. 
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APPROXIMATIONS 

The complete simulation of the reactor in three dimensions is 
clearly impracticable, since there are about 4000 channels each 
containing about 7 fuel elements. The principal points of interest 
as regards reactor safety and output are the axial temperature and 
flux distributions in the channels. A steady-state solution of the 
equation which describes an imaginary homogeneous reactor 
yields a method of averaging the radial variations of flux and 
shows that the reactor output can be described in terms of an 
average channel. We therefore simplify the problem to the con¬ 
sideration of the simulation of one channel, with reactor flux a 
function only of axial position in that channel. 

The fuel elements introduce axial variations into the exact 
solution for the flux for a homogeneous reactor, which is a cosine 
function. However, these variations can safely be averaged, and 
are small compared to some other approximations which must be 
made if a simulation of reasonable size is to be possible. The flux 
equation is thus a partial differential equation in axial distance z 
and in time t. Similarly, the irregularities in the temperature dis¬ 
tribution in the radial direction in the channel, moving outwards 
from the centre of a fuel element, are neglected. For the simula¬ 
tion to be possible at all, the partial differential equations for heat 
transfer and for flux must be approximated by finite differences. 

A further uncertainty is the behaviour of the graphite modera¬ 
tor. Both its thermal conductivity and its specific heat vary with 
time in the reactor in ways which are not fully understood, and 
must therefore be guessed on the basis of existing knowledge. 

However, in spite of these and other approximations, which can 
all be justified and the errors resulting from them estimated, an 
extremely useful simulation can be set up. 

DESCRIPTION OF THE EQUATIONS 

The equation which describes the neutron behaviour, the neu¬ 
tron flux equation, is given as equation (1). It is a neutron 
balance equation saying that the rate of loss of neutrons by axial 
flow (in an element of length dz) plus the rate of gain of neutrons 
from the reactivity and the delayed neutrons gives the rate of 
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change of the number of neutrons. The reactivity term is com¬ 
posed of the natural reactivity, the reactivity control, and the 
reactivity effects of temperature. The delayed neutrons are 
emitted after a certain half-life by fuel atoms which have absorbed 
a neutron a short time before. 

The flux equation is a partial differential equation in both space 
and time. Since we are interested primarily in the time variations 
of the reactor, we split the equation into a set of finite difference 
equations in space. For a simulation of reasonable size, that is, 
below about 200 amplifiers, the flux equation is written in finite 
difference form as four equations giving the flux at four planes in 
the reactor. 

At these four planes, heat-transfer equations describing radial 
flow of heat from the fuel and graphite to the coolant are written 
down. The axial flow of heat can be neglected, since the tempera¬ 
ture gradients in the radial direction are very much greater than 
those in the axial direction. The flow is taken as symmetrical in 
the angular direction. 

The heat produced by the nuclear reaction is proportional to the 
neutron flux. The greater part of this heat is produced in the fuel, 
but a small proportion appears in the graphite, can and gas. The 
heat-transfer equations (3) to (7) give heat balance of the form: 

Heat generated by flux, minus heat conducted to the can, 
plus heat conducted from the centre fuel, per second, 
equals rate of change of heat in the fuel surface. 

The coolant equations (8), (9) give the heat balance for the 
channel divided into four parts centred on the four planes of the 
simulation. These equations are of the form: 

Heat leaving the section, minus heat entering from the pre¬ 
ceding section in the gas, equals heat gained from the fuel 
can and graphite. 


The Flux Equations 
The basic equation 


be written in the form 

+ {fc oo (l-jS)-l}.^ + 5k oo 0 

+k (—¥- + —-—+ l 0 D<l>, 

+ 1+D/rJ 9 Zt 09 
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where <j> = neutron flux, z measures distance along the channel, 
k x = reactivity constant of reactor. 

$k 0 o = reactivity control, temperature and xenon effects, 

P = /?'+/?", Ml, l 0 are constants. 


K\ 

D 


P' ! P" \ , _ delayed neutron flux 
1 + D/A' 1 + D/A'y' 9 _ representation, 

_ d . pS' . , 

= —, and is the neutron source term. 


We now write the term in d 2 (f>/dz 2 in the form of a set of four 
finite difference equations. First the length of the reactor core is 
divided into eight equal lengths L/2 numbered as shown in Fig. 31 
from 1 to 9. The flux at any plane i is then written as </>;. 

If the three-point finite difference approximation in the form 
4>i = 4>i+i + <f>i- 1 ~2</> f is used, errors of up to 10 per cent can 
arise. Five-point Lagrange expansions about points 2, 4, 6 and 
8 are therefore taken. In order to satisfy the physical boundary 
conditions on the value of <£, at planes 1 and 9, planes 0 and 10 
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are introduced, where the flux falls to zero. These planes are a 
distance a 2 L below plane 0 and a 0 L above plane 8. Further, two 
planes -2 and -8, with “imaginary” fluxes of -<p 2 and -</> 8 , 
are taken symmetrically below plane 0 and above plane 10. Then 
the flux in the neighbourhood of plane 2 is written in terms of 
-tj ) 2 , cj ) o = 0, (j ) 2 , <j ) 4 and $ 8 ; and near plane 4, in terms of 
(j) 0 = 0, (j> 2 , 0 4 , (j ) 6 and (j> s ; and similarly for planes 6 and 8. 
Using this method, it can be shown that 


2 /W = 21 + 12a 2 -4a 2 4a 2 (3-a 2 ) 

\dz 2 ) 2 2(1 + 3a 2 + 2a 2 ) l + 3a 2 +2a 2 

a 2 (3 2a 2 ) 


2(2 + 3a 2 + a|) 

+ 3+2-2 3 + Oj 

\cz 2 / 4 3a 2 l + a 2 2 + a 2 


3(3+ a 2 )- 


w . 

- 1 A, 4 

3 + a 8 


th 


3(3 + a 8 )'^ 4 

2 + a 8 

1 + a 8 
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1 + 3 a, 




+ 

3a 8 

_ 

a 8 (3-2a 8 ) 


4a 8 (3-a 8 ) 

rh 

v^ 2 A" 

2(2 + 3a 8 + a| 

j • ^4 + 

1 + 3a 8 +2a§ 

• 


21 + 12a 8 —4a 8 
2(l + 3a 8 + 2a|) ^ 

Exercise 

Derive these expressions using the methods given in Chapter 


Heat Transfer Equations 
For the fuel rod and can: 

(vpc) fl DT ft = (hA) fl/2 (7> 2 - T ri ) + S fl <j>, (3) 

(vpc) f2 DT f2 = (hA) flfl (T f -T f2 ) 

-(hA) f ' S (T f2 -T s ) + S f2 <t>. (4) 

For the fuel can to coolant: 

(vpc) s DT s = (hA) fiS (T fl - T s ) - (hA) s JT s - T e ) - 1 f K. (5) 
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For the graphite: 


(vpc) mi DT mi = ( hA) mimi (T mi - T mi ) 

- (hA) muC (T mi — T c ) + S mi <t> + l c R, (6) 

(vpc) mi DT mi = (hA) mimi (T m ~ T mi ) + S mi <f>, (7) 

where 

vpc = volume x density x specific heat (cal), 
hA = heat transfer coefficient x area (cal/°C sec), 

S = heat produced per unit of flux (cal/unit flux sec), 
l c R = heat radiation from can to graphite (cal/sec), 


and where the suffixes f u f 2 , s, c, m u m 2 refer to the regions of 
fuel can, coolant, and moderator shown in Fig. 32, and a pair of 
these suffixes refers to the heat transfer between the two regions, 
across their common interface. 



Coolant Equations 

Consider the heat balance of the gas in the region of plane 2, 
that is, between planes 1 and 3. Using the notation of the previous 
section, 

(hA\' C (T s - T C2 ) + (M) m>c (T mi - T C2 ) 

= W.c p .(T C3 - T C1 ) + D(M. c p . T C2 ), (8) 
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where T , T C1 , T C3 are the coolant temperatures at planes 1, 2 
and 3, W is the mass flow in g/sec, c p the specific heat of the gas, 
and the mass of gas in the region is 


M = 


mass flow 
gas velocity 


length of region = WIJV. 


Assuming that the coolant behaves as a perfect gas, it is easily 
seen that the term M.c p . T Ci is independent of temperature, since 
the temperature and density are inversely proportional. 

We assume a linear rise of temperature from T Cl to T C3 ; thus 
the equation can be written as 


, (MU 


T _(hA) s (hA) m , 


T ml + T cl . 


This form is convenient, since the heat transfer depends on 
mass flow of gas through the channel, and can be written in the 
form 

where h 0 refers to the heat transfer coefficient at full power mass 
flow of W 0 , and f(W/W 0 ) is an experimental relation, approxima¬ 
ting to a square root law. 


Exercises 

1. By writing the coolant equation in the form indicated, show that 
one function generator and a multiplier are needed in the simulation 
of the coolant. Draw out the simulation circuit for this equation, and 
for the equation in the form of (8) above, where the heat-transfer 
coefficients are taken to depend on mass flow as given later in the 
paragraph. What accuracy troubles are to be expected at low mass 
flow for these two methods of simulation? What cost is paid in 
amplifiers and equipment for a four-plane simulation to obtain the 
I improved accuracy? (A simulation circuit is shown in Fig. 33.) 
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2. If Ts and T mi ' are the absolute temperatures in degrees Kelvin 
of the can and moderator surface, the radiation term R is 

R = 2nR [in . ff(7V 4 — T m ,' 4 ) (cal/sec unit length), 

where R [in is about 2-5cm, and a is Stefan’s constant = 1-32 x 10 12 
c.g.s. units. If T s and T mi lie between 200° and 600°C, give circuits and 
scaling factors for generating /?, = (Ts ' 4 — T,, H ' 4 ) using 

(a) function generators, 

(b) by expanding as a power series in Ts and r m „ 

(c) by factorizing, generating the factors and then multiplying up. 

What are the faults and how accurate can each circuit be expected 
to be? 

(Hint-. T,' = Ts + t; take (T.l 5), (Tmi 5) and (3/J/1010) as a first 
try.) 



Fig. 33a. 






Rack 1 


Potr. No. 

Gain 

Input 

Output 

AO, A8, 

0 00508 

~Hi 

15 

A1 



A2 




A3 




A4 




AS 




A6 

002768 

15 

15 

A7 

002768 

19 

15 

A8 

See AO 



A9 




BO 

0-2200 

02 

07 

B1 

2-305 

13 

12 

B2 

0-1229 

19 

12 

B3 

0-8771 

(-^) 

12 

B4 

0-2200 

12 

07 

B5 

1-1984 

09 

10 

B6 

0-5992 

10 

10 

B7 

0-0032 

R 

10 

B8 



B9 

! 



CO and C9 

Not used 

Initial condition 


DO 




D1 

0-1137 

-Hi 

08 

D2 

0-3225 

09 

08 

D3 

0-3225 

08 

08 

D4 

0-4612 

10 

09 

D5 




D6 

0-1336 

+Hi 

09 

D7 

0-2652 

08 

09 

D8 

1-188 

09 

09 

D9 



EO 




El 

0-00877 

19 

19 

E2 

0-00877 

15 

19 

E3 

0-00065 

R 

19 

E4 



E5 




E6 

0-03507 

06 

19 

Rack 2 

AO-9 

Not used 



BO 

1-0 

20 

03 

B1 




B2 




B3 

2-563 

20 

05 

B4 

2-563 

05 

05 

B5 



B6 

0-2348 

20 

06 

B7 

0-2348 

06 

06 

B8 

0-03946 

05 

20 

B9 

0-01336 

06 

20 
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GENERAL PROCEDURE FOR PREPARING 
A SIMULATION 

SIMULATION CIRCUITS 

Once the basic forms of the equations have been fixed, and all 
reasonable simplifications and approximations made, sketch 
circuits should be drawn out. These must be simplified as much 
as possible, and a preliminary estimate of the amount of equip¬ 
ment required, together with any special circuits can be given. 

The circuit for one plane of the example is given in Fig. 33. 

PHYSICAL DATA 

The next job in the preparation of the problem is the gathering 
of accurate data for all the equations. In the example, the thermal 
capacities of the fuel rods, moderator and can must be found. 
Design figures for the size of the core, the behaviour of the control 
rods, and the reactivity at each plane must be dug out of the 
engineers and physicists responsible. Decisions on which sets of 
data are to be used to obtain a self-consistent simulation, and 
how the uncertainties of the data are best handled, must be made. 
This is usually a lengthy and exacting job. 

As soon as a rough set of values has been found, the equations 
can be scaled and the rough gains given. The normal procedure 
of guessing at maximum values and re-scaling will give finalized 
circuit diagrams, except for the precise values of gains. 

For a large programme of work on the simulation, it is well 
worth while writing or having written a simple digital computer 
programme of the interpretive type to give the gains for the 
simulation from the physical data. A more complex programme 
will give the expected steady state and other information. As soon 
as accurate data are available, the machine gains and steady states 
can quickly be found; and when later experiments give altered 
thermal conductivity, heat-transfer rates and so on, the computer 
saves a great mass of tedious calculation. The digital programme 
also removes the worry of faulty arithmetic, and its printout can 
give the simulation gains in a neatly tabulated typescript form 
rather than a possibly untidy handwritten sheet of paper. 
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PATCHING AND SETTING UP 

Once the equations have been scaled and the circuits fixed, the 
problem can be patched. Any special circuits should be checked 
on the machine or on a test computer—e.g. the simulation of the 
radiation from can to moderator might be checked for errors and 
noise. 

The main layout of the problem on the patchboard should be 
worked out before patching starts, so that leads are kept reason¬ 
ably short. As the patching proceeds, check off each amplifier, 
multiplier or potentiometer on the circuit diagram and equipment 
lists. All subsequent alterations and additions should be entered 
carefully on these lists, and the lists and circuits kept with the 
tabulated values of gains. As each section or subsection of the 
problem is wired on to the patchboard, check it through com¬ 
pletely to make sure that the wires go to the right places. Cases 
have been known of amplifier outputs connected to earth, 
potentiometer inputs and sliders inverted, multiplier potentio¬ 
meters connected upside-down and so on, in spite of the efforts 
of the patchboard designer to make things easy for the user. 


PROCEDURE ON THE MACHINE 

When the patchboard is fully wired up and all amplifier lists, 
equipment lists and circuit diagrams are ready, the problem should 
be set up on to the machine systematically, marking each gain as 
it is set. Machine faults will make themselves evident slowly— 
some part of the circuit will be found to integrate rapidly off to 
saturation due to a broken lead, or a multiplier will refuse to 
track round properly. It is as well to recheck the patchboard, 
gains and amplifier switching before calling angrily for the main¬ 
tenance man. The operator who can say, “The slider of potentio¬ 
meter B6 seems to have fused” and it has, is more popular with 
the maintenance men than the one who says, “I can’t get anything 
from potentiometer B6—see to it” and he has forgotten to patch 
any volts to it. Accurate and careful work on the machine saves 
time for the operator, the maintenance staff, and most important 
saves the cost of machine time. 
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It should be possible to obtain some sort of steady state for new 
problems involving 200 amplifiers in a fortnight or so of hard 
work. 

Once the operator has a steady state which he trusts, the pro¬ 
gramme of runs can be started. As an extension to the steady 
state, a standard run can be chosen and used each day to check 
that the machine is functioning correctly. A machine can run 
perfectly for weeks, and then suddenly go haywire for no apparent 
reason. Eventually it will be found that, for example, a capacitor 
unit has been changed in the course of routine maintenance, and 
it has been incorrectly set up, or it has an internal fault. 


PROGRAMME OF WORK 

The programme of runs required, sets of different gains for 
parameter surveys, and development work required should be 
prepared as the simulation is set up. Runs can then be taken from 
the machine at the highest rate, since everyone on the machine 
knows exactly what is required. Active production time, as 
distinct from preparation, will be divided about equally between 
preparing for the run, and checking up on the machine immedi¬ 
ately prior to the run, and then taking the run. In between batches 
of production, gains will have to be reset and circuits changed, and 
a certain amount of machine idle time is almost inevitable as the 
runs from the machine are investigated. 


SOME TYPICAL SIMULATION RESULTS 

Two typical runs from a reactor simulation are shown in Figs. 
34 and 35, re-plotted for clarity. If one of the four gas circulating 
blowers fails for some reason, the gas mass flow through the 
reactor will fall rapidly from 100 to about 75 per cent, whilst the 
reactor remains at constant power. The curves in Fig. 34 show 
the sector control rod movement to reduce that power, and the 
rise in T sm , the maximum fuel can temperature. T sm does not rise 
above about 480°C, so the sector rods move fast enough to 
prevent a dangerously high fuel temperature being reached. The 
third curve shows the variation in gas outlet temperature from the 
reactor. 
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The second set of curves, in Fig. 35, show the reactor behaviour 
on loss of all blower power due to voltage failure to the circulator 
motors. T sm rises to show 515°C, and this is sufficient to trip the 
reactor, and causes all control rods to fall into the reactor, thus 
rapidly reducing power. The axial temperature profiles inter¬ 
polated from the results at different times after the faults are shown 
in the figure. 


EXAMPLES OF THE USE OF LARGE 
ANALOGUE COMPUTERS 

THE NUCLEAR POWER INDUSTRY 

An example from the nuclear industry is described in some 
detail earlier in this chapter. Simulations such as the one described 
are widely used in safety studies throughout the nuclear industry. 

Other uses to which analogue techniques have been applied are 
studies of the build-up and decay of the fission product xenon, 
which has a very marked effect on the operation of the reactor as 
it is a strong absorber of neutrons. In a reactor operating at 
power, the element iodine is created as a fission product at a 
constant rate. The iodine then decays to xenon with a half-life 
of about 8 hr, and the xenon decays and is destroyed by the 
reactor flux to form caesium. At power, the xenon reaches a 
steady equilibrium level, but if the power is reduced to zero, the 
xenon is no longer destroyed by the flux but is still being created 
by decay from the iodine, and the concentration of xenon there¬ 
fore increases. Xenon is a strong absorber of neutrons, and so the 
reactivity of the reactor is greatly reduced for a period of about 
16 hr, and a return to power is often impossible. Many aspects 
of this “xenon poisoning” problem have been simulated, from 
simple simulations of the average xenon concentration following 
a load change to complex studies involving up to 300 amplifiers 
to investigate the radial stability of the flux distribution in the 
reactor. It is possible for a local disturbance in the flux to build 
up in any horizontal plane into waves of increasing amplitude, 
due to the time lag between the xenon and the flux. One method 
of solution is to write the equations in modal form, like the 
different modes of vibration of a stretched string. The equations 
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to describe each mode can then be solved, and the flux obtained 
by summing the modes. Such studies have shown that the 
reactors being built in Great Britain for power stations must be 
controlled as a set of independent subsections, each automatically 
regulated by its own control rods. 

A large amount of work has been done on the non-nuclear 
aspects of nuclear power stations, such as the heat exchangers, 
control rod drives and associated control loops, the electrical 
stability of the transformers, switching gear, motors and gener¬ 
ators under normal and fault conditions and the simulation of a 
complete nuclear station in designing and developing automatic 
control loops for the control of the power station considered as a 
single unit. 

In short, the nuclear industry offers a wide range of problems 
which are susceptible to analogue means of solution. 


AIRCRAFT SIMULATION 

The aircraft and missile industry makes a great deal of use of 
analogue computation. Apart from the extensive use of com¬ 
puters in the investigation of the automatic controls required to 
move control surfaces, vary fuel flow, steer radar and direction 
finding aerials, the autopilot system and many others, the actual 
flight dynamics are frequently simulated. 

By choice of coordinate system, the equations of motion of an 
aircraft can be written down fairly simply from Newton’s laws of 
motion. But the equations are seriously non-linear and complex— 
there are six degrees of freedom of an aircraft, three of velocity 
and three of rotation. 

The equations are frequently linearized for small disturbances 
about the steady level-flight condition. The reader is probably 
familiar with the technique of replacing a variable x by x 0 + Sx, 
where x 0 is the steady-state value. All terms in the second or 
higher orders of small quantities are neglected as small, the 
steady-state terms cancel, and a set of equations are left which are 
linear, with the deviations from the steady state as variables. 

A further simplification is to ignore the coupling between the 
rool, yaw and pitch motions about x-, y- and z-axes, taking the 
x-axis along the aircraft. Experiment shows that this coupling is 
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often slight. Studies based on these simplifying assumptions are 
often valuable as an aid in the design of the aircraft and its 
associated control systems, but many refinements can be added to 
the basic 5 or 6 degrees of freedom simulation. The modes can be 
coupled and the complete large angle displacement equations, 
involving rotation of axes transformations in the form of direction 
cosine matrixes or Euler angle relationships can be simulated. 
The equations for the drag terms can be refined and the automatic 
systems and auto-pilot behaviour added, until highly complex 
systems result. 

Stability studies and parameter surveys are two uses for such 
simulations. Once an aircraft is built, a simulator may be set 
up to match its response, and then used for training pilots, using 
a complete cockpit and instrumentation display system tied to the 
computer simulation. 


MISSILE SYSTEMS 

The behaviour of a guided missile system has been the subject 
of a great deal of analogue investigation, much of it secret. The 
systems for data gathering, radar searching, tracking, ground- or 
missile-based guidance or control, together with the airframe 
dynamics and the rocket motor behaviour have all been simulated, 
both individually and together as systems. The advantages of such 
simulations are typical of the advantages of the analogue approach. 
The behaviour of a costly, complex, ill-understood system can be 
evaluated prior to its construction. Optimum values of the para¬ 
meters can be evaluated, and parameter surveys can be under¬ 
taken. Many of the aspects of the interrelation of the many linked 
servo-controls and control loops can only be evaluated on an 
analogue machine. Feasibility studies can be done for future 
systems. 

A particular problem in missile systems is the response to noisy 
input information. The adjoint methods (Chapter 6) can be used 
to reduce the large number of computer runs needed to evaluate 
the response of a time varying system to noise. The optimum 
firing conditions for the minimum target “miss” distance are 
required under many conditions, and the adjoint method gives as 
output the root mean square miss distance in a single run. 

8 ACM ' 
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A technique which has been used is to fix on a particular com¬ 
puter run of importance, and to devote considerable effort to 
obtaining a digital computer programme to give the system 
behaviour for that run. The comparison of the analogue simula¬ 
tion output and the digital answer should give added confidence 
in the results of the many other types of run which the simulation 
can produce. The established digital computer answer will serve 
as a check whenever any alterations or additions are made to the 
analogue system. 


PROCESS CONTROL 

The investigation of chemical plant can be very complex due to 
the large number of variables involved. In practice, design is very 
reliant on empirical methods, experience and extrapolation from 
previous results. Even approximate mathematical expressions for 
the dynamic behaviour of the plant are very hard to obtain. 

Analogue computers have been used with success in many 
cases, and their application is increasing. If possible, the equations 
describing the process are written down and simulated in the 
ordinary way. But the amount of equipment can quickly become 
prohibitive. Simplifications must be made by linearizing equa¬ 
tions, taking a coarse set of finite difference points for partial 
differential equations, ignoring small terms and so on. 

An approach which has been used for a distillation process is to 
write down the large number of equations describing the evapora¬ 
tion and condensation on about a hundred separate plates in the 
distillation column and to evaluate the transfer functions relating 
each input to each output in turn on a digital computer. Assum¬ 
ing that the process is linear for a range about the design condi¬ 
tions, the system can be simulated by simulating these transfer 
functions. The resulting simulation may only involve three or 
four amplifiers for each transfer function, and it can then be used 
to study the best design of the control loops to maintain the out¬ 
put by varying the flows and temperatures in the system. One of 
the problems peculiar to chemical plant simulation is the accurate 
production of exponential functions which govern the reaction 
rates. Diode function generators can be used, but other circuits 
have been developed for the various special aspects of this prob- 
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lem. One technique is to transform all the equations into logarith¬ 
mic form. This can be used to simulate systems in which a very 
wide range of variables must be scaled. Various automatic re¬ 
scaling devices have also been tried with some success. 

An important difficulty in simulating process plant is the 
adequate simulation of a distributed system such as a bank of 
heat transfer tubes in a heat exchanger. The plant parameters 
can be averaged, and the whole process assumed to take place at 
one point, or finite difference models can be set up. But an 
accurate representation requires a large amount of equipment, 
and can be difficult to manage on the computer. 

Associated with such distributed systems is the difficulty of the 
transport delay. The length of time taken by a fluid to flow from 
one end of a pipe to the other cannot be simulated properly. The 
transfer function is of the form 


e-'^l -,T + ^ 3 


p 3 t 3 


The approximations are given in Fig. 24 (p. 68). 

Various other methods have been tried in simulating the trans¬ 
port delay, using tape recorders with movable heads, scanning 
uniselector systems, ordinary strip chart recorders with photo¬ 
electric readers and so on. The results are seldom very good, and 
they are difficult to use. 

In general, the advantages of the analogue computer in chemical 
process plant studies are those indicated previously. The standard 
techniques are normally sufficient for most problems, but large 
computers are often needed. The increasing use of analogue 
methods promises design and construction economies, and 
releases experimental plant for other work. 


CONTROL ENGINEERING 

Many of the techniques used in simulation on electronic ana¬ 
logue computers come direct from control engineering. Indeed, 
in some ways electronic analogue computers have been developed 
because of the power of the analogue in the field of control 
investigations. Many problems involve non-linearities and com¬ 
plexities far beyond the capacity of paper and pencil analysis 
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using the most powerful tools which the studies of control theory 
have given the specialist. The very nature of the analogue lends 
itself directly to the simulation of transfer functions—the gain of 
an amplifier with input and feedback impedances Z L and Z 0 is 
G(P ) = -Zq/Zi, and complicated impedances can be built up out 
of resistors and capacitors to give a wide range of transfer 
functions from input to output of the amplifier. (Chapter 1, Fig. 
2a (p. 13), gives a short list of useful impedances.) In addition, 
non-linear processes can often be simulated with ease using 
multipliers and diode circuitry. 

The advantages of speed of solution, ease of parameter altera¬ 
tion, ease of programming and low cost have been mentioned 
elsewhere in the book. Actual parts of the system involved can 
sometimes be incorporated for testing into an analogue investiga¬ 
tion, e.g. a control rod and its motor can be included into a small 
simulation of a nuclear reactor. Simulations are often used in 
optimization and in feasibility studies of control systems for such 
reasons. The stability and accuracy of control of a system by its 
automatic or human operator plus automatic control can be 
determined in advance of its construction by the design engineer. 
And simulations of control systems can be used to train operating 
personnel. Link trainers for aircraft and the Calder Hall reactor 
simulator are good examples of this use. 

A technique which is of wide use in control studies is that of 
repetitive computation. The simulation is time scaled to allow 
solution rates of 50 or 60c/s, and relay or other switching arranged 
to give a cycle of set initial conditions—compute—clamp—set 
initial conditions, at the 50 c/s rate. The solution is then displayed 
on a cathode ray tube, and, if required, photographed for record 
purposes. On such a simulation, the effect of altering the system 
parameters—loop gain, or the amount of phase advance, for 
example—can quickly be seen. This mode of solution is often 
less accurate than the conventional “real time” simulation, but 
has the advantage of speed and simplicity of solution. 

The fields of application of analogue computers to control 
engineering are too great for more than a very few examples to 
be mentioned. Simulation has been used in steel rolling mill 
design, for the servos which work the rolls; in the design of 
feedwater controllers for boilers and reactor plant; in the control 
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and design of surge tanks for water turbines; in the design of 
aircraft engine controls and of the servos to move the aircraft 
rudder and positional controls; and in practically every field 
where anything is controlled. 
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chapter 6 


GENERAL 

ANALOGUE 

DEVICES 


THE DIFFERENTIAL ANALYSER 

The principle of the planimeter was extended by Prof. J. 
Thomson, the brother of Lord Kelvin. He published a paper in 
the Proceedings of the Royal Society in 1876 “On an integrating 
machine having a new kinematic principle”. He noted that the 
major fault of the planimeter was the sliding action of the wheel 
and, working on an idea of Prof. J. Clark Maxwell, produced the 
device shown diagrammatically in Fig. 36. A hardened disc can 
rotate about its axis, which is inclined at a slight angle to the 
vertical. A cylinder is mounted so that it can rotate about its 
axis, which is horizontal and perpendicular in direction to the axis 
of the disc. A ball rolls on the disc, resting against the cylinder, 
the whole being arranged so that the point of contact of the 
ball and disc lies on the horizontal line through the centre of 
the disc. 

Suppose now that the ball, by means of a cage arrangement, is 
held with its centre a distance y from the centre of the disc. If the 
disc rotates through a small angle dx from its angular position x, 
the ball will transmit a rotation of y dx to the cylinder. The total 
rotation of the cylinder as x and y vary will then be proportional 
to Jjdx. 

In the same issue of the Proceedings of the Royal Society, Prof. 
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Fig. 36. 


Sir W. Thomson, later Lord Kelvin, shows how to use this 
instrument to solve a second-order differential equation 

d_/l du\ _ 
dx\Pdx)~ U ’ 

which, quoting from a Cambridge examination question of 1874, 
is “easily solved” thus: 

Try as a first guess the function u = u y as solution. 

Form ii 2 = J*P^C-j'*ii 1 dj^dje 

u 3 =J P^C-J u 2 dx^dx, etc. 

The series of functions u u u 2 , u 3 ,... converge to a solution of 
the given differential equation. 

Now Prof. Thomson formed these integrals with the aid of his 
brother’s machine until he found a solution. “But then came a 
pleasant surprise. Compel agreement between the function fed 
into the double machine and that given out by it...” by establish¬ 
ing a connection from output to input (in modern terms a feedback 
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connection), and “The motion of each (integrator) will thus 
necessarily be a solution. Thus I was led to a conclusion which 
was quite unexpected; and it seems to me quite remarkable that 
the general differential equation of the second order with variable 
coefficients may be rigorously, continuously, and in a single 
process solved by a machine.” 

From this pleasing conclusion Lord Kelvin proceeds in a second 
paper to show how to solve the general differential equation with 
variable coefficients using the integrating device. He envisaged its 
use for the solution of many problems involving vibration such as 
a non-uniform stretched cord, a wave in a non-uniform canal, the 
conduction of heat and the differential equation of the tides. 

This device was not developed. However, in the 1920’s V. Bush, 
working at M.I.T., independently developed an integrating device 
consisting of a wheel resting against a disc. A small hardened 
wheel with its axis horizontal rolls on a horizontal disc. The disc 
is arranged to be able to rotate about its axis and also to move 
on a carriage beneath the wheel in the direction of the axis of the 
wheel. As before, if the centre of the disc is positioned a distance y 
from the plane of the wheel by the carriage, and the disc rotates 
through an angle dx, the wheel rotates through an angle ydx, and 
so the shaft of the wheel gives an output of j>dx as x and y 
vary. 

Using this principle, a six-integrator machine was built, accurate 
to about OT per cent, the connection between the shafts being 
made via gear trains. 

This machine was followed during and after the war by several 
other machines, ranging from a two-integrator machine built out 
of Meccano by Prof. Hartree, and found to be accurate to better 
than 1-5 per cent, to a thirty-integrator machine at M.I.T., 
accurate to one part in 10,000 with paper tape controlled set-up 
and initial conditions for input, and accurate photoelectric 
couplings between the shafts, connected by relay circuits. 

Problems such as the behaviour of emitting filaments, the 
control and stability of synchronous motors, and propagation 
along transmission lines were solved on these machines. 

The torque available on the integrator output shaft is limited by 
the friction of the wheel against the disc; it must therefore be 
amplified in some way. This can be electrical, but has usually 
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Ouptut arm 



been achieved by mechanical torque amplifiers using counter¬ 
rotating drums with friction bands controlled by the input to 
give power to the output shaft. A study of Fig. 37 will indicate 
the action of this device. 


SET-UP OF A DIFFERENTIAL ANALYSER 

The integrator has been described; it is represented in Fig. 
38a. The input to the disc position is w, and v is its rotation, 
giving u = fvvdu as output. Inputs to the machine are provided 
from an input table by manually directing a pointer over a pre¬ 
drawn graph of the required input q =f(p), a function of the 
variable p. Thus q might be a forcing function input to a system 
as a function of time, taken as p. Summation and multiplication 
by a constant is accomplished in gear boxes, shown in Fig. 38c. 
And an output table is shown in Fig. 38d; a pencil is moved 
across a paper according to the relationship between y and x 
given by the machine. A simple set-up gives the method of 
solution of an equation of the type 

d 2 y dy 

£! + /■£ + «)’-H*). 

This is shown in Fig. 39. 

Note that multiplication can be performed by two integrators 
together, using the relation uv = $udv + $vdu. Layouts may be 
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derived for the generation of many functions using the property 
termed “regenerative connection”, or, in electronic analogue 
terms, implicit function generation. However, the method is of 
much wider application, as integration can be performed with 



L 

o*- 

q P 
q = ftp) 

(b) 



p/q = n/m 


(el 

Right hand Left hond 
(e) 



Fig. 38. (a) Integration unit, (b) Input table, (c) Sum¬ 
mation and constant multiplication, (d) Output table, 
(e) Gearing between longitudinal and cross-shafts. 


respect to any variable with the mechanical integrator, whereas 
the electronic methods of integration rely on the charging of a 
capacitor with time, and therefore give integration with respect 
to time only. This method has found an application in incremental 
computers, also called digital differential analysers, which are 
briefly described later. 
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Exercises 

1. Derive layouts to generate the functions 
d 2 y 

sin 6 and cos 0 from — + y = 0. 

1 r \ 1 

log e x and - from J - dx = log* x, - d(log* x) 


2x 


2. Show that any gear system with an infinitely variable ratio can be 
used as an integrator. 


THE INCREMENTAL COMPUTER 

One of the advantages of the differential analyser is its ability 
to perform integration with respect to any variable in the com¬ 
puter. For some applications this facility with very great long-term 
stability are essential, for example in the construction of flight path 
and guidance-computing equipment for aeroplanes and missiles. 
Various aircraft firms in America and in England have built equip¬ 
ment for this purpose which has been referred to under the names 
“incremental computer” and “Digital Differential Analyser” or 
DDA. 

Unlike a true analogue, the variables are represented by binary 
digits and their increments. A diagram is shown in Fig. 40a. A 
variable y is stored in a register Y, to which may be added its 
increments d y by an adding unit. Increments dx in a variable x 
open a gate which allows y to be added to a register R, called the 
remainder register. Thus y is added into I? at a rate depending 
on the rate at which increments in x arrive, so that R accumulates 
a number proportional to the integral of y . dx. Now R will fill 
up progressively and will overflow. The overflow or spill digit 
from R can thus be used as an output ds-k.y.dx. Various 
different refinements to this basic scheme allow both signs of input 
to be accepted by, for example, taking an input of zero to be a 
subtraction and of 1 to be an addition, so that in fact a steady level 
in the value of a variable is represented by a train of l’s and 0’s. 
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Improvements of technique and in the mathematical understand¬ 
ing of the process has enabled very high accuracy integrating units 
of great stability to be built. These can then be interconnected to 
solve differential equations, almost as if they were analogue com¬ 
puting units, in association with adding and other equipment. Two 
conventions for drawing out the integrating unit are shown in Fig. 


(o) 



Increments in 
ydx spill from R 

fcicrements in x 


Increments in y 





Fig. 40. (a) Block diagram of DDA element, (b) Conventional 
representations of DDA element, (c) Two simple DDA circuits. 


40b. Since integration with respect to variables other than time is 
possible, the units can be interconnected in exactly similar circuits 
to the layouts of differential analyser equipment to give, for ex¬ 
ample, multiplication by using the fact that d(u.v) = u.dv+v.du, 
and other implicit functions as explained for the differential 
analyser. 

A point of importance is the relation between speed of solution 
of a problem set-up and the accuracy of solution. For a given type 
of circuitry, the speed of the adding and digit transfers will be 
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fixed. But, if the variable y is represented by a ten-binary digit 
register, the accuracy will be half that of an eleven-binary digit 
representation; however, the speed of solution (determined by the 
fineness or coarseness of the increments dy in the least significant 
digit and the fixed rate of transfer of digits) will be twice as great 
for the ten-digit case. 

With modern techniques, incremental computing equipment can 
be built which give solution speeds competitive with conventional 
analogue equipment, with very superior accuracy, but at about 
twice the cost. Transistorized equipment gives great reliability and 
small size, and this, together with the high accuracy, have been 
the features particularly required in aircraft systems rather than 
low cost. Incremental computing equipment has no great use 
elsewhere in the analogue field, due to its high cost. 

Exercise 

Repeat the examples given for the differential analyser in incremental 
computer terms. 


CONDUCTIVE ANALOGUES 

The solution of Laplace’s equation V 2 <f) = 0 and of Poisson’s 
equation W 2 4> = 4np in one form or another are basic problems 
in engineering. For complicated boundaries, the exact solution 
becomes impractical. One of the techniques which has been widely 
used in their solution is the analogue of the flow of an electric 
current in a conducting medium. 


resistance sheets 

A large number of engineering problems can be reduced to two- 
dimensional terms or to problems of cylindrical symmetry. For a 
two-dimensional field satisfying the Laplace equation, a uniform 
thin sheet of copper was used in 1845 by Kirchhoff. However, the 
conductivity of copper is too great for accuracy, and a convenient 
substance was not produced until the Western Union Telegraph 
Co. introduced “Teledeltos” paper in 1948. This paper is made 
by adding carbon black to the paper pulp during manufacture and 
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then coating one side with a thin laquer, the other with a very thin 
layer of aluminium paint. The resulting paper is about 0 004 in. 
thick with a resistivity of the order of 20 kQ measured from edge 
to edge of a square of paper of any size. Its resistivity varies by 
about 10 per cent in any sheet, and it is different for different 
directions of current flow, but it is cheap and easy to use. Bound¬ 
aries may be marked by using a conducting silver paint, and it can 
be cut to shape with scissors and stuck to a backing plate for 
strength. A blunt metal sensing probe, or a ball-pointed probe 
can be used to measure the field pattern resulting from the bound¬ 
ary voltages by balancing the probe voltage in a bridge circuit. 
An accuracy of about 10 per cent is to be expected, which can 
often be within the limits of the problem data. The speed and ease 
of solution, even with such limited accuracy, made this a useful 
engineering tool. 

For a solution of adequate accuracy, it is necessary that the 
conductive sheet is uniform in resistivity and is isotropic. Clearly 
local variations in resistance, systematic variations and different 
values in different directions will give false results, and it is these 
variations that are hardest to overcome satisfactorily. The elec¬ 
trodes used for the boundaries must have a low resistance com¬ 
pared to the sheet, and the input impedance of the probe must be 
high compared to the resistivity of the sheet or measurement in¬ 
accuracies will result. Another factor is the coefficient of resistivity 
with temperature, which must be low to prevent change of resist¬ 
ance due to the heating of the sheet from the current passing 
through it. The humidity of the atmosphere can also have a 
marked effect on the resistivity of the sheet. 

Specially made conductive sheets have been used as well as the 
Teledeltos paper. One method is to spray a layer of conductive 
paint on to paper, another is to weave a cloth from alternate 
strands of metal and fibre. Rubber loaded with carbon can be 
used. One fairly accurate method is to use ordinary transparent 
draughtsman’s paper, which appears to have no fibre structure 
and a very high degree of uniformity. Its conductivity appears to 
be due to absorbed water and salts, which form a uniform electro¬ 
lyte throughout the material. The resistivity of the paper is about 
3 x 10 10 f2 per square, it varies widely with humidity, and it makes 
very carefully designed sensing equipment necessary. But con- 
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ducting boundaries can be drawn on the paper with an ordinary 
pencil. 

For the solution of Poisson’s equation, the continuous source 
term must be approximated by discrete current sources at various 
points on the paper, fed with appropriate voltages from resistance 
chains. Another method is to use the property of superposition 
of the fields, and use only one source electrode, moved from point 
to point, and the fields resulting from each point are then added 
together. 

The sensing probe can be very simple or highly elaborate, with 
lead screws to position the probe, mounted on a carriage above 
the paper. The probe must make a good contact with the paper 
without damaging it and without wearing its point flat on the 
paper by repeated movements across the surface of the sheet. A 
simple Wheatstone bridge arrangement with a meter to indicate 
the null can be used to plot out the equipotentials point by point. 
Alternatively, valve voltmeter devices or electrometer valves for 
high resistance sheets can be used to obtain the position of the 
field lines. 


THE ELECTROLYTIC TANK 

The electrolytic tank is a method of great power in the solution 
of Laplace’s equation. The first such system is described by W. C. 
Adams in 1875, and it has been widely used since then. As in the 
conducting sheet analogue, a model of the system is made by 
enclosing in a tank of electrolyte a scale model of the boundaries 
of the problem. Equipotentials are made of metal, and insulating 
material is used for the streamlines. The model is excited by a.c. 
at between about 300 c/s and 10 kc/s and a probe dips into the 
electrolyte to measure the voltage. 

It has been found that ordinary tap-water is often adequate as 
an electrolyte, though a low concentration of copper sulphate or of 
acid in distilled water has some slight advantages. As with the 
conductive sheet, the resistivity must be uniform, linear and of the 
right order of magnitude compared to the electrodes and probe. 
The electrolyte must be clean, it must not react chemically with 
the electrodes, and its surface tension should be low to minimize 
the effect of the meniscus at the probe tip. An a.c. excitation is 
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chosen of a frequency which represents a compromise between 
the polarizing effect of d.c. at the electrodes and the increasing 
capacitive effects at higher frequencies. 

By using a square wave excitation, both these effects can be 
effectively eliminated. Typically, about 20V excitation between 
electrodes is convenient. 

For problems in two dimensions, simulating the behaviour of a 
medium of constant properties, a simple tank with the equi- 
potential boundaries simulated by appropriate conductors will 
give the required field patterns. However, if the properties of the 
medium vary, for example the varying dielectric constants of the 
ceramic, paper and air insulation of a transformer bush, bound¬ 
aries are necessary within the tank. One method is to change the 
depth of the tank in a step, but vertical conducting wires are 
needed at the step to make the current densities approaching and 
leaving the step uniform in the depth of the tank, whilst following 
the correct field lines in the horizontal plane. Another method is 
to represent the different media by separating different electrolytes 
by thin barriers with a large number of vertical wires to transmit 
the current from one to the other on the lines of the above. 
Ordinary blotting paper separators and porous sponges in the 
electrolyte to increase the resistivity are other methods which have 
been used. 

For problems with cylindrical symmetry, a tilted tank can be 
used. A wedge-shaped quantity of electrolyte now forms the 
model, and the above methods of simulating changing physical 
properties can be used. 

The simplest type of probe is a metal wire just dipped into the 
liquid. But this is subject to errors due to the meniscus of the 
liquid against the wire. Multiple arrays of probes can be used 
close together, and the resulting signals used to steer the probe 
array along the equipotentials. Various servo systems have been 
designed to plot the field lines automatically, but the complexity 
needed is not usually warranted, since the greater part of the 
problem time is taken in setting up and checking the problem. 
However, they reduce errors and tedium, and changes of the 
electrolyte with time have less effect on the results since the read¬ 
ings are taken much quicker. 

With a very carefully designed system, and using great care in 
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all measurements, accuracies of 0-1 per cent have been obtained. 
However, 1-2 per cent is an average figure to be expected from 
normal equipment. 

The method has been widely used for the design of electron lens 
systems, for insulation problems in high voltage equipment, and 
in various problems of heat conduction. Some experimental work 
has been done to determine its possible usefulness for three- 
dimensional problems, but without wide application. 


RESISTANCE MESHES 

A method for the solution of field problems that has developed 
rapidly since the war is the resistance mesh. The continuous 
medium of the electrolytic tank or conductive sheet analogue is 
replaced by discrete resistances in a rectangular mesh and cur¬ 
rents are fed in at the mesh boundary or at interior points. The 
voltages at the mesh nodes are read off using high impedance volt¬ 
meters, and the equipotentials can then be plotted out. 

We have seen that we may approximate to the solution of the 
equation 

d 2 f , d 2 f n 


by taking four equally spaced points 1, 2, 3, 4 around a fifth point 
number 0, and defining/ 0 from/i,/ 2 ,/ 3 , and/ 4 by 

4/o=/l+/ 2 +/3+/4- 


Suppose now that four points 1, 2, 3, 4, are joined to the point 0 
by equal resistors R, and the voltages V u V 2 , V 3 and V 4 are pro¬ 
portional to f u f 2 , fi and / 4 . Then we write down the current 
balance for the node 0 as 

R R R R 


Thus 4 Fq = V 1 + V 2 +V 3 + V 4 . 

Thus, if the area of interest in the field problem is divided up 
into a rectangular mesh, and a mesh of resistors built up to link 
the nodes of this mesh together, when the boundary values of the 
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problem are replaced by proportional voltages at the boundaries 
of the resistance mesh, the node voltages give a finite difference 
analogue to the problem. The procedure is as follows. First 
visualize a conductive sheet analogue of the original field, made 
of a material of resistivity p per square, divided up by a rect¬ 
angular mesh of spacing 8x and 8y, shown in Fig. 41a. 



Fig. 41. 


The resistance mesh analogue is shown in Fig. 41b. In this, the 
resistor carries current in the ^-direction corresponding to the 
current which would flow in the sheet analogue in the y-di recti on 
between AB and CD. 

So the resistance of R A must be proportional to the resistance 
between AB and CD, which is p(8y/8x). Similarly, in the x- 
direction, R 2 is made proportional to the resistance between PQ 
and RS, that is, p(8x/5y). Thus, we associate each resistor with 
a “vector area” of the field. In an exactly similar way, in three 
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dimensions a vector volume is associated with each resistor of a 
three-dimensional mesh. If the mesh spacing is 5x, 8y, 8z, this 
gives resistors 

_ 8x 5y Sz 

x ~ p JiTz' Ry ~ p d^c’ Rz = p bTsi‘ 

At boundaries the resistors can be calculated in exactly the same 
way. For a square mesh at a boundary parallel to the mesh the 
vector area is half that of the other elements, so that the resistor 
must be twice the value of the others. Values for other boundaries 
are shown in Fig. 42. 
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It is possible to use a finer mesh to cover areas of special 
interest by matching the boundaries and deriving the resistor 
values by the area rule as before. 


Exercises 

1. Calculate resistors for a resistance mesh analogue of the current 
flow across a square sheet ABCD from the side AB at voltage V to 
CD which is earthed, with a circular hole in the centre. If possible, set 
up the mesh and test it. 
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2. Set up a resistance mesh to study the edge effects of a parallel 
plate condenser. 

3. Derive the resistor values for a problem in cylindrical polar co¬ 
ordinates using the vector volume method. Set up an analogue for the 
steady heat flow in a piston head from the flat face to a hemispherical 
cooling cavity, assuming no heat loss from the cylinder walls. 


EXTENSIONS OF THE SIMPLE RESISTANCE MESH 


A great deal of ingenuity has been used in devising methods for 
the solution of other types of problem to those described by the 
Laplace equation. For example, Poisson’s equation can be simu¬ 
lated by feeding currents to the nodes from external sources. This 
can make the setting up lengthy as each time a source current is 
fixed, due to the resistance of the various external resistors which 
give the voltages for the boundary conditions, everything alters. 

Dr. G. Liebmann and others have published details of various 
special resistance meshes for the solution of equations such as 
V 2 (j>+k4> = 0 (in which the input currents to the nodes are 
adjusted on an iterative procedure towards the solution), for the 
solution of the diffusion equation V 2 0 = kd<j>ldt (by expanding 
as a finite difference equation in both space and time) and for the 
solution of the wave equation V 2 (j> = k(d 2 4>jdt 2 ). 

One of the most ingenious is the double mesh for the solution 
of problems in elasticity. The mathematical expressions defining 
problems in elasticity are hard to formulate and harder to solve, 
and the correct specification of boundary conditions can be very 
complex. The basic equation is the biharmonic equation 


, 8 4 <j) d A 4> a 4 4> n 

V dx 4+ 8x 2 dy 2 + dy 2 ~°- 


This may be solved by cascading two resistance meshes one 
above the other, the upper one solving for an auxiliary stress 


function 


8 2 4> 8 2 (j> 
8x 2 8y 2 


which is specified along the boundaries of the problem and the 
lower solving for $ which can be the displacement. The nodes of 
the upper mesh are connected to those of the lower by resistors; 
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it can be shown that if these resistors are large compared to the 
mesh resistors, the mesh voltages form the required analogue. If 
resistors feed in currents to the upper mesh nodes, an analogue 
of V 4 (j) =/, which also arises in elasticity, is produced. Satisfac¬ 
tory results have been obtained from such networks using mesh 
resistors of about 100Q, interconnecting resistors of 100-200 kQ, 
with the two meshes working at about 10V and 2V. 

Dr. Redshaw and Dr. Rushton of Birmingham University have 
developed a negative resistance, using a simple transistor unit, to 
simulate various other problems in elasticity involving buckling 
and flexure. For the deflection ca under a tensile load P and a 
lateral load q, the differential equation involved is of the form 


d 4 co 

d^' 


Writing 
this becomes 


d 2 co 

M = -El 


In finite difference form, this becomes 


and 


h 2 

£»1 + ©3 2C00 — 0 . 


These can clearly be solved by a resistance mesh for a tensile 
load, but a negative resistance is needed for a compression to 
provide for the term (Ph 2 /EI)M 0 , which becomes negative. 
Various problems have been solved using these negative resistance 
units. They are cheap, compact and easy to set up, and provide 
a method of overcoming the costly and cumbersome nature of 
inductance-capacitance networks which can also give negative 
impedance network elements. 

Practical use of the resistance mesh has been made in many 
engineering problems in elasticity such as the bending of beams, 
the torsion of bars and hollow sections and the bending and 
extension of plates. Other uses are the many examples in engineer¬ 
ing where the Laplace equation occurs, for example the flux 
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distribution in electrical machines, electromagnetic fields in cavity 
resonators, forces on transformer windings and the steady state 
neutron flow in a nuclear reactor. 


RESISTANCE-REACTANCE AND OTHER MESHES 

In a precisely similar way to the set-up for resistance meshes, 
capacitors or inductors can be included in meshes to simulate 
other types of system. For example, the heat conduction equation 
is of the form V 2 (f> = k(d<j>jdt). If each node of the resistance mesh 
simulating the space derivatives in V 2 </> is connected to earth by 
a capacitor, we simulate immediately the thermal capacity of the 
system and thus the whole differential equation. In one dimension 
this analogue is illustrated by the analogue between the equation 

0i + 02-2<Ao = 

at 


ViV 2 


ZVo_ r dV 0 
R dt ' 


The reader may readily see how this may be extended to two or 
three dimensions, and how the rules for finding the values for 
resistors are extended for the capacitor values. By suitable time 
scaling, the simulation can be arranged for repetitive display on 
a cathode ray tube, or for real time computing. Values of node 
voltages can be recorded using high input impedance recording 
devices. 

In a similar way, using the fact that the current in an inductor L 
is (l/L)j(F 1 —F 0 )d/, we may simulate the wave equation using 
inductance-capacitance meshes, and, by adding resistors, the 
modified wave equation 


, d 2 d> dd> 


The understanding of the details of analogues of these two 
equations is left as an exercise, as it follows immediately from the 
above. 

There are a large number of mixed inductance, resistor and 
capacitor mesh analogues which have been devised, but few have 
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been used with any great success, due to the high cost, bulk, and 
poor accuracy of such devices in practice. In elasticity, various 
transformer mesh systems have been devised, but these are for 
the specialist. The leakage effects when using inductors, together 
with their internal resistance and stray capacity, limit the accuracy 
from practical meshes to about 5 per cent, although the mathe¬ 
matical theory developed by G. Kron in the course of his studies 
of analogues can be very powerful when a digital machine is 
available. 


SOME SPECIAL APPLICATIONS 

MONTE CARLO METHODS 

The mathematical analysis of a process subject to statistical 
laws can be very complex in some cases. A problem of some im¬ 
portance in nuclear engineering is the determination of the fine 
structure of neutron flux in a reactor fuel element and channel. 
Digital methods are normally used, the machine storing the 
information as to the channel geometry and neutron behaviour. 
The machine, typically, will consider a neutron of a certain 
velocity, allow it to move for a random length of time, make a 
collision and bounce off at a random angle and repeat with 
another random path, until a third random number indicates that 
the neutron is absorbed at a collision. The numbers of paths, 
collisions and their locations are stored and then printed out after 
sufficient have been computed to provide a statistically valid 
result. 

An analogue computer has been built to simulate this process 
by the English Electric Co. For a neutron of velocity V, moving 
in a channel with X- and T-axes of symmetry at an angle 8 to 
the X-axis, its position is determined by 

x = J Fcos0dt + x o , 

y=J Fsin0df + y o - 

Thus x and y can be generated by integrators from V and 8. 
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Since the geometry of the reactor can be built up from repeated 
reflections in the boundaries of the unit cell (Fig. 43), a neutron 
passing from one channel to the next may be simulated by 
reflecting it at the boundary of the cell. This may be accomplished 
by detecting when the X- or T-coordinate reaches the boundary 
value, and then changing the sign but not the magnitude of the 
relevant component of velocity to reverse the direction of integra¬ 
tion. Fast transistor switching circuits are employed for this 



purpose, in association with standard computing amplifiers 
developed for the computer Saturn described on p. 80. 

The random time between collisions is generated from a small 
radioactive source with a scintillation counter. On receipt of a 
pulse from the counter, the integration is stopped and an angle 
generator started. The angle generator produces sinO and cost?, 
and is stopped after a time determined by a pulse from another 
scintillation counter, giving the direction of the reflected path. 
The process then repeats, until an absorption is indicated at some 
collision by another device. The region in which the absorption 
occurs is determined, and the absorption is counted automatically 
by decatron counters. 

Many refinements can be added to simulate different properties 
of different regions and to detect the positions of the absorptions. 
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An optical technique for simulating asymmetrical configurations 
for shielding and other calculations has been perfected, displaying 
the neutron position on a cathode ray tube and detecting its 
position by photoelectric means. It has been found that 20-30 
paths can be computed per second, which is comparable to the 
fastest digital computers using advanced programmes. 

It is possible to use the computer to compute the complete 
neutron behaviour from its very fast initial velocity when born in 
a fission, through the complex slowing down process, to the 
thermal velocities at which it can cause fission. Using the com¬ 
puter in this way, shielding problems involving very complex 
geometry can be solved to a high accuracy. Due to the direct 
optical method of representing the geometry, problems impossible 
to solve on the fastest digital machines can be readily handled. 


THE ADJOINT METHOD 

Mathematically it may be shown that the response of a linear 
system to a noise input may be determined by a knowledge of the 
noise characteristics together with the system frequency response, 
or of the response of the complete system to an impulse. For a 
system which does not change with time, it is usually more 
laborious to find the frequency response experimentally than the 
impulse response. 

However, for a system which changes with time (e.g. a rocket 
whose mass decreases with time), the response to a noise input 
depends both on the instant of observation t 2 and also on the 
time t 2 -t u during which the noise input has existed. Mathe¬ 
matically this makes both the frequency and the impulse response 
methods lengthy and laborious. 

A technique which has been evolved to simplify the procedure 
makes use of a transformation of the independent variable of the 
differential equation which describes the system. The resulting 
adjoint set of equations can give a simulation which shows the 
system behaviour in a reversed time direction. Simulation time 
starts from the moment of observation and runs backwards; the 
simulation gives the system response at time t 2 to a noise input 
existing from t x to t 2 , in terms of f, as independent variable, the 
time on the simulation. 
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The mathematics of this process is beyond the scope of this 
book, and is described in Roger and Connelly, Analog Computation 
in Engineering Design, McGraw-Hill. The following is a brief 
description of the method. 

If the impulse response of the system is h_ 1 (t 2 ,t 1 ), for observ¬ 
ation at time t 2 following an impulse at time t u then h \(t 2 t,) can 
be non-zero only where t 2 > t x . That is, in the shaded region of 
lug. 44a. If, then, an impulse is applied at t u the history of the 



Fig. 44. Response following paths AB and CD. 


system is given by following a path such as AB, as t 2 increases 
from zero (Fig. 44b). The adjoint system follows paths such as 
CD, m the direction of the arrow, where x = t 2 -t l is simulation 
time. The output which is normally required from the system is 
the root mean square response to noise. This is obtained by 
squaring and integrating the ordinate of the path CD with respect 
to T. 

Forming the Adjoint Circuit of a Simulation 

Having found the adjoint simulation of a given simulation, the 
following rules have been shown to apply by Lanning and Battin. 

1. Turn each element of the circuit round, and reverse the 
direction of signal flow. 

2. Make all time varying parts of the system start from t 2 and 
run backwards in time. 

3. Interchange outputs and inputs of the system. The new 
input will be an impulse <5 (t) at time t, and the output will 
be the impulse response /i_, (/ 2 , — T ). 
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4. Integrators and summers have gains of 1 only, and potentio¬ 
meters have one output only. The gains may be reappor¬ 
tioned after the transformation. 

An example given by Rogers and Connelly is: 

1 d 2 0 o 1 dfl 0 
Kf(t) d t 2 + KM dt o_0i 

the equation giving the output 0 o of a time varying system to 
input 9 t . 

The simulation circuit is shown in Fig. 45a. 

Applying the above rules, the adjoint simulation is shown in 
Fig. 45b. 



Fig. 45. 


In general, any term in the original equation of the form 
f„(t)(d n /dt n ) is replaced by a term 

and the time t of observation is replaced by t 2 - 
Thus the equation given above becomes 

I il . 1 ! 

Xdr 2 [ f(t 2 - t) J Kdr{ f(t 2 -r ) J 

+ h_i(/ 2 ,f 2 —t) = (5(r). 

This is the equation solved by the adjoint circuit above. 
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OPERATIONAL RESEARCH 

A great deal has been written about the use of analogue 
computers in operational research, but it seems that little practical 
work of importance has been done to date. 

Basically, the type of problem which is solved is the determin¬ 
ation of the optimum conditions of operation of a system, from 
a knowledge of its restraints expressed as inequalities. A diode 
simulates an inequality directly and summing amplifiers will 
simulate the linear relations between the variables of a linear 
problem. 

As an example, it is desired to minimize (or maximize) a 
variable y (perhaps the cost of a product) where y is a linear 
function of variables x 1 ,x 2 ,x 3 ,...,x n (costs of the operations 
required to produce the product). 

Thus y = a 1 x 1 +a 2 x 2 + ... +a„x„, 

where some or all x, ^ 0 (i= 1,2,...,«). 

Further, there are constraints on the variables of the form 
A t = fliiX 1 +a 12 x 2 + ...-l-a lB x B gR 1 
A 2 = o 2 i x l + a 22 x 2 +... + a 2n x 2 S R 2 

A m = a ml X l+ a ml X 2 + ■ ■ -+ a mn X n = Rm 

(For example, the limits on the costs of research, development, 
etc.) 

For constant values of oq and a u this type of problem is called 
a linear programming problem. 

We use a method of steepest descent—that is in this case, find 
the gradient of the function y as a function of the x ; and move 
along that direction until a constraint is reached. 



where i 1 ,i 2 , — ,i„ are unit vectors in the directions of the x 1; x 2 , 
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Thus, starting from any point we must subtract 

(or add) increments a. t 8, a 2 0, a n 9 from the x-coordinates in 

order to optimize y. In addition, these increments must be 
constant with 

A l ^ R u A 2 ^R 2 ,...,A n <:R n . 

We therefore fix a velocity for each x f by 
dx* 

—± = 0a * — K £ dji e p summing from j = 1 to j = m, 

where e } - Aj-Rj if A s > Rj, 

= 0 if Aj <■ Rj. 

The X; will therefore move freely until y “runs into” a boundary 
defined by some R Jt then move a little beyond the boundary until 

a t 9 = Kjj'Ojtej, 

and then continue along the boundary. The point therefore 
moves until y is optimized. In order to minimize y, 6 is taken 
negative; to maximize, 9 is taken positive. The relative values of 
8 and K determine the rate at which the solution is approached 
and the accuracy of solution. For the simulation, ej are generated 
from summers, with a diode feedback to give zero output if 
Rj ^ Aj. 

The values of Ke 1 ,Ke 2 ,...,Ke m and the a f 0, taken as inputs to 
integrators, then give x,. A diode limiter prevents the values of 
X; becoming negative. 

This type of technique can give solutions to problems of the 
nature—what is the best way of allocating labour to a job, or of 
setting up shops to cover a marketing area. One problem of 
importance is the allocation of loads to power stations on an 
electrical supply grid to minimize the distribution costs and to 
ensure that the efficient stations are used in preference to the 
inefficient and costly stations. It is possible to handle non-linear 
problems in a limited way, but the danger in all optimizing meth¬ 
ods is that a false minimum is found before another “deeper” 
minimum is seen. 

Related types of problem are the analysis of the rise in sales of 
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a product following an advertising campaign, or of the personnel 
requirements of a firm with fluctuating sale and staff recruitment 
rates. It is likely that the development of these analyses will effect 
considerable economies in some fields. 
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The field of analogue computation ■ 
and simulation is covered in this : 
book in a compact and handy form 
Those reading the book are assumed 
to have a knowledge of the calculus, 
mathematics and physics such as 
would suffice to gam admission to a 
university. A certain amount ol 
mathematics is presented and the 
components and construction of an 
analogue computer are discussed 
: Simple examples aregiven of its use, 
followed .by- a detailed discussion of 
a larger problem. Finally, other : 
i techniques, such as conductive 
analogues, are discussed, and some 
■Typical applications of electronic 
analogue computers are given. The 
book aims at avoiding excessive 
detarl while yet giving the prospective 
; user a concise grounding in the use cf 
: an :angl.Qgue.c'omputer, Typical uses 
of analogue machines are described, 
in -addition to the basic principles. 
Examples are given wherever possible, 
Technical colleges.that give courses 
in analogue methods will welcome 
this book and it will be. invaluable . . 
for engineers and management 
wishing-to gain an appreciation of 
the scope of analogue computation 






